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1.0  INTRODUCTION  AND  SUMMARY 


1.1  BACKGROUND 

Future  DOD  large  space  system  concepts  such  as  High  Altitude  Large 
Optics  (HALO),  High  Energy  Laser  Optics  (HELO)  and  Millimeter  (MM)  Wave 
applications  have  stringent  1 ine-of-sight  (LOS)  pointing  and  jitter  perfor¬ 
mance  requirements  which  cannot  readily  be  met  with  existing  technology. 

Under  the  Active  Control  of  Space  Structures  (ACOSS)  Program,  the  Defense 
Advanced  Research  Projects  Agency  (DARPA)  has  sponsored  a  number  of  studies 
for  the  development  of  a  unified  structural  dynamics  and  control  technology 
base  to  support  the  future  development  of  large  space  systems  missions. 

Major  emphasis  In  these  studies  has  been  placed  on  generic  control  law  develop¬ 
ment  for  active  vibration  suppression. 

This  document  reports  the  results  of  one  of  these  studies,  namely 
ACOSS-8,  conducted  by  TRW  Defense  and  Space  Systems  Group  for  DARPA  from 
May  1980  through  April  1981.  The  study  addressed  three  areas  of  active 
structure  control/vibration  control: 

•  Stability  ensuring  designs  in  the  presence  of  modal  truncation/ 
inaccurate  structural  models  using  a  method  based  or.  the  posi¬ 
tivity  of  operators. 

•  On-orbit  parameter  identification  using  the  maximum  likelihood 
method. 

•  Tradeoffs  between  model  reference  adaptive  coitrol  and  fixed 
gain  control  with  and  without  on-orbit  parameter  identification 
and  system  tuning. 

1.2  SUMMARY 

One  of  the  challenges  in  the  active  control  of  large  space  structures 
(LSS)  is  to  assure  stability  of  control  system  designs  based  on  inaccurate 
structural  models.  The  two  main  errors  ir.  the  dynamic  model  of  a  structure 
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arise  from  modal  truncation  and  inaccurate  knowledge  of  the  modal  frequencies 
and  mode  shapes  obtained  by  finite  element  methods.  Additional  errors  are 
introduced  by  poor  knowledge  of  structural  damping  and  the  effects  of 
structural  nonlinearities. 

What  is  needed,  therefore,  is  a  stability  ensuring  design  methodology 
that  does  not  rely  on  modal  truncation,  has  only  low  sensitivity  to  the  exact 
knowledge  of  the  structural  model  (robust),  and,  with  some  modifications,  can 
be  extended  to  nonlinear  systems.  We  are  presenting  such  a  method  here, 
based  on  the  positivity  of  operators.  This  method  yields  good  performance, 
but  it  is  not  necessarily  optimal  when  significant  model  uncertainties  are 
present,  because  robustness  with  respect  to  knowledge  of  the  modal  data 
tends  to  occur  at  the  expense  of  a  more  conservative  design. 

A  major  result  of  this  study  ir.  a  design  philosophy  for  active 
control  of  LSS  that  can  provide  robust  stability  and  also  meet  stringent 
performance  requirements.  In  order  to  assure  robust  stability  the  positi¬ 
vity  method  is  applied  to  design  the  active  structural  controllers.  Once 
the  entire  system  is  on-orbit,  tests  will  be  conducted  to  identify  critical 
model  parameters  (e.g.  modal  data  and  structural  damping).  Parameter  identi¬ 
fication  methods  which  appear  most  suitable  for  this  are  the  maximum  likeli¬ 
hood  method  and  some  of  the  frequency  response  methods.  In  the  present 
study  we  concentrated  on  the  maximum  likelihood  method.  Data  frcm  the  on- 
orbit  identification  tes  ,s  a.:  transmitted  to  ground  where  tne  identification 
algorithm  is  applied  and  the  a  priori  knowledge  of  model  paramaters  updated. 
The  improved  model  parameters  are  then  used  to  recompute,  on  the  ground,  new 
control  parameters  in  order  to  tune  the  system  to  achieve  the  desireo  perfor¬ 
mance.  Figure  1-1  illustrates  the  concept. 

An  alternative  to  explicit  parameter  identification  and  subsequent 
system  tuning  is  direct  adaptive  control.  Model  reference  adaptive  control 
(MRAC)  is  particularly  attractive,  since  it  allows  specification  of  the 
desired  performance  through  a  reference  model.  Figure  1-2  shows  a  basic 
model  reference  adaptive  configuration.  Among  the  advantages  of  MRAC  are 
(1 )  autonomy,  (2)  automatic  on-line  controller  adjustments  for  component 
failures,  and  (3)  automatic  on-line  controller  adjustments  in  response  to 
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Figure  1-1.  On-Orbit  Model  Identification  and  System  Tuning 


changing  conditions,  e.g.  changes  in  the  thermal  environment  or  reorienta¬ 
tion  of  large,  articulated  sub-structures.  On-line  implementation  of  MRAC 
is,  however,  more  complex  than  parameter  identification  and  subsequent 
system  tuning.  Tradeoffs  between  these  control  approaches  have  been 
conducted  in  this  study. 
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Figure  1-2.  Dioital  Model  Reference  Adaptive  Control 
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STUDY  CONCLUSIONS 


The  positivity  approach  has  been  shown  to  work  well  for  the  design 
of  LSS  control  Systems.  It  has  been  successfully  applied  to  the  design  of 
the  attitude  control  and  active  structural  damping  systems  for  the  Draper 
Model  #2  whose  payload  LOS  was  inertially  stabilized.  The  resulting  control 
system  performed  well  and  was  robust  with  respect  to  model  parameter  varia¬ 
tions.  To  achieve  optimal  performance  as  specified,  the  system  had  to  be 
subsequently  tuned,  assuming  improved  knowledge  of  the  system  model  (mode 
shapes,  frequencies  and  structural  damping).  In  the  positivity  approa*'*'  the 
improved  model  knowledge  allows  then  the  use  of  more  elaborate  "embeddin-j 
operators"  (related  to  the  classical  Popov  multipl iers  and  absolute  stability 
sectors)  which  leads  to  less  conservative  stability  results,  e.g.  higher 
gain/bandwidth,  and  therefore  better  performance.  If  the  system  can  justi¬ 
fiably  be  assumed  linear,  the  tuning  process  can  directly  apply  the  multi- 
variable  frequency  domain  design  techniques  of  MacFarlane  and  others,  using 
frequency  domain  shaping  of  the  open  loop  characteristic  gains  of  the 
system. 

As  the  study  results  have  clearly  shown,  there  exist  tradeoffs 
between  the  robustness  of  a  design  and  its  performance:  one  is  obtained  at 
the  expense  of  the  other.  Model  uncertainty  may  be  large,  but  to  generate 
a  robust  design,  the  accuracy  of  the  LSS  model  must  be  known  within  some 
specifiable  limits.  If  the  model  uncertainty  is  small,  robustness  may  be 
traded  for  performance. 

Model  uncertainty  is  reduced  through  on-orbit  parameter  identifica¬ 
tion.  Using  the  maximum  likelihood  method,  excellent  estimates  of  physical 
and  modal  LbS  parameters  were  obtained.  Structural  damping  is  the  most 
difficult  parameter  to  estimate  because  it  is  generally  very  small  so  that 
the  measurement  data  exhibits  only  a  low  sensitivity  to  the  parameter  and 
longer  data  lengths  are  required.  Furthermore,  the  mathematical  parameter¬ 
ization  of  structural  damping  is  usually  not  accurate,  i.e.,  not  truly 
representative  of  the  physical  phenomenon.  The  estimation  of  physical  LSS 
parameters,  such  as  mass  distribution,  cross-section  of  members,  dimensions, 
polar  moments,  spring  constants,  etc.,  instead  of  modal  parameters  has 
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significant  advantages,  since  it  drastically  reduces  the  number  of  para¬ 
meters  that  must  be  estimated. 

Model  reference  adaptive  control  (MB  )  is  an  alternative  to 
explicit  on-orbit  parameter  identification  „.id  control  system  tuning.  In 
this  study  MRAC  was  applied  to  control  a  perturbed  model  of  the  Draper 
tetrahedron  structure  (Draper  Model  #1)  and  very  good  results  were  achieved. 
MRAC  was  shown  to  be  a  viable  approach  to  LSS  control  and  it  works  well. 
However,  the  study  showed  that  model  parameter  identification  followed  by 
a  controller  tune-up  results  in  still  better  performance,  because  the 
human  designer  can  use  his  more  global  knowledge  of  the  system  in  the  design. 
For  example,  poor  observability  of  critical  modes  ..i  the  sensor  measurements 
can  be  accounted  for  (to  a  certain  degree)  by  a  designer,  while  a  model 
reference  adaptive  controller  is  completely  oblivious  to  the  situation. 

This  was  the  case  in  the  tetrahedral  structure,  where  the  modes  which  are 
critically  affecting  the  motion  of  the  apex,  are  only  poorly  observable/ 
controllable  from  the  supporting  bl-peds  of  the  tetrahedron  where  the 
sensors  and  actuators  are  located. 

The  conclusion  is  that  explicit  on-orbit  model  parameter  identifica¬ 
tion  and  subsequent  controller  tuning  via  ground  command  is  the  preferred 
approach.  Model  reference  adaptive  control  should  only  be  used  in  cases  .where 
spacecraft  autonomy  is  absolutely  essential  and/or  where  unexpected,  non- 
preprogrammable  changes  in  the  system  are  likely  to  occur.  Otherwise,  the 
increased  on-board  software  complexity  of  model  reference  adaptive  control 
cannot  be  justified. 

Finally,  we  conclude  that  the  proposed  three-step  design  approach 
for  LSS  control  systems  appears  to  be  the  correct  approach  in  the  sense 
that  it  minimizes  risk  with  the  potential  for  maximizing  performance. 

1)  It  assures  stability  wich  robustness  via  the  positivity  design 
method. 

?)  It  identifies  model  parameters  through  on-orbit  tests. 
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3)  It  tunes  the  controller  to  maximize  performance. 

1.4  RECOMMENDATIONS  FOR  FURTHER  STUDY 

Furcher  studies  in  the  positivity  design  approach  should  address 
example  designs  with  non-collocated  actuators  and  sensors,  nonl inearities 
in  the  controller,  and  further  work  on  the  development  of  an  easy-to-follow, 
computer-aided  design  methodology. 

In  on-orbit  parameter  identification,  non-zero  initial  dynamic 
conditions  of  the  spacecraft  should  be  addressed.  Frequency  domain  identi¬ 
fication  techniques  should  be  studied  in  more  depth  and  compared  with  the 
maximum  likelihood  method. 

Finally,  for  realistic  design  applications,  controller  implementa¬ 
tion  issues  must  be  addressed.  Foremost  among  these  are  actuator  and 
sensor  development,  modelling,  and  laboratory  test. 

1.5  ORGANIZATION  OF  REPORT 

Section  2  addresses  the  positivity  design  approach  and  Section  3 
on-orbit  parameter  identification.  Model  reference  adaptive  control  and 
performance  comparisons  with  fixed  gain  controllers  are  presented  in 
Section  4.  Quoted  references  are  listed  following  Section  4.  Appendix  A 
contains  a  description  of  the  Draper  Model  =12. 


sj 
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2.0  THE  STABILITY  ENSURING  DESIGN  APPROACH  VIA  POSITIVITY  \ 

The  stability  ensuring  design  methodology  makes  use  of  positivity/  | 

energy  dissipation  theory  to  ensure  the  stability  of  the  LSS.  The  key  results  i 

are: 

•  The  negative  feedback  connection  of  a  positive  system  with  a 
strictly  positive  system  is  stable 

•  An  LSS  with  collocated,  ideal  (by  ideal  we  mean  infinite  band- 
with),  actuators  and  rate  sensors  is  a  positive  system 


•  Embedding  operations  permit  the  treatment  of  non-positive 
systems,  e.g.,  systems  with  realistic,  bandwidth-limited 
actuators  and  sensors. 

The  usefulness  of  the  theory  becomes  apparent  as  one  realizes  that  no 
assumptions  are  made  about  the  dimensionality  or  linearity  of  the  systems 
involved.  It  is  for  this  reason  that  the  theory  permits  the  treatment  of 
modal  truncation,  stability  robustness  and  nonlinearities.  The  price  one 
pays  for  the  generality  is  that  designs  tend  to  be  conservative. 

It  has  therefore  been  proposed  to  use  positivity  in  conjunction  with  on- 
orbit  identification  and  control  "tune-up".  Positivity  will  be  used  in  the 
pre-flight  design  where  there  is  parameter  uncertainty;  then  as  information 
is  gathered  on-orbit,  multivariable  frequency  domain  techniques  combined 
with  positivity  (to  treat  non-linearities)  can  be  used  to  tune  the  design 
for  high  performance. 

It  is  clear  that  as  system  know^dge  increases,  the  need  for  a  conser¬ 
vative  pre-flight  design  diminishes.  In  cases  where  good  pre-flight  informa¬ 
tion  is  available  over  approximately  twice  the  control  band,  the  positivity,/ 
multivariable  frequency  domain  approach  may  be  applied  directly  without  on-orbU 
identification. 

The  positivity  method  has  been  successfully  applied  in  this  study  to 
various  examples:  a  flat  plate  (Section  2.2)  the  Draper  tetrahedron  (Section 
4.2.2)  and  the  Draper  Model  =2  (Section  2.3).  In  each  case,  the  positivity 
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approach  led  to  a  very  simple*  and  stability  robust  design  with  good  perfor¬ 
mance.  Subsequent  controller  tuning  based  on  assumed  improved  plant  modelling 
(obtained  via  on-orbit  parameter  identification,  as  in  Section  3)  then  per¬ 
mitted  all  designs  to  meet  all  the  desired  performance  specifications. 

2.1  The  Positivity  Design  Methodology 

2.1.1  General  Background  On  Positivity 

The  positivity  method  is  a  direct  outgrowth  of  V.  M.  Popov's  work  on 
absolute  and  hyperstabil ity.  The  basic  underlying  theory  is  general,  makes 
extensive  use  of  functional  analysis,  and  is  not restricted  to  linear  time- 
invariant  plants.  In  fact,  the  method  has  in  the  past  been  associated  with 
the  stability  of  nonlinear  systems  (Lur'e  problem)  where  both  the  nonlinear 
and  linear  parts  of  the  system  must  meet  certain  positivity  conditions. 

These  conditions  are  only  sufficient  conditions  for  stability  and  therefore 
tend  to  be  conservative  at  times.  Since  well-known  necessary  and  sufficient 
conditions  exist  to  establish  the  stability  of  linear  systems,  it  is 
apparent  why  positivity  techniques  have  in  the  past  generally  not  been  applied 
to  linear  systems. 

The  stability  problem  in  the  control  of  large  space  structures  (LSS) 
differs  from  lore  conventional  linear  systems,  because  the  plant  to  be 
controlled  is  theoretically  of  infinite  order  (and  practically  of  very  high 
order),  while  the  control  1 er/estimator  must,  of  necessity,  be  of  relatively 
low  order.  It  appears,  therefore,  desirable  to  use  tests  which  do  not  require 
explicit  evali  ;t ion  of  the  closed  loop  system,  but  which  impose  conditions 
individually  on  the  plant  and  the  controller,  assurino  stability  when  the 
loop  is  closed  by  negative  feedback.  The  stability  theory  based  on  the 
positiyity  of  operators  is  a  powerful  method  having  just  this  characteristic. 


*  Typical  frequency  domain  compensation  consisted  of  second  or  fourth  order 
lead-lag  filtering. 
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Positivity  Theorem:  The  system  S  (in  Figure  2.1)*  is  asympto¬ 

tically  stable  if  at  least  one  of  the  loop  transfer  matrices  G  and  H  is 
strictly  positive  real  and  the  other  is  positive  real. 

The  main  advantage  of  this  result  is  that  conditions  are  imposed 
individually  on  the  plant  G(s)  and  the  controller  H(s),  ensuring  stability 
when  the  loop  is  closed  by  negative  feedback. 

For  ease  of  understanding,  and  because  they  are  presently  of  main 
interest,  we  have  considered  here  only  the  linear  stability  result.  However, 
the  proof  f  the  Theorem  given  in  Reference  [1],  treats  the  controller  and 
the  plant  as  general  operators,  not  necessarily  linear,  defined  on  a  Hilbert 
function  space.  This  extension  of  the  theory  to  nonlinear  systems  will  be 
summarized  and  used  in  Section  4. 


VECTOR  VALUED  SYSTEM  "S" 

Figure  2.1  System  Considered  By  The  Positivity  Theorem 


2.1.2  Definition  And  Meaning  Of  Pc:itive  Real 

Positive  real  means  exactly  what  it  says.  The  real  part  of  a  transfer 
function  as  function  of  frequency  is  always  greater-than-or-equal-to  zero. 
Strict  positive  real  means  the  real  part  is  always  greater  than  zero. 

For  single-input-single-output  (SISO)  transfer  functions  g(jio), 
(strict)  positive  realness  means  therefore  that 

Re  |  g(joj)|  >  (>)  0,  for  all  ~<x>  <(^<*0 

For  multiple-input-multiple-output  (MJMO)  systems,  the 
*  The  theorem  is  also  true  when  H  is  in  the  feedback  path. 
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definition  needs  to  be  more  precise  as  one  is  now  dealing  with  a  transfer 
matrix.  Here,  a  square  MIMO  transfer  matrix  G(j<u)  is  (strictly)  positive 
real  i f 


G(jw)  +  G*(jw)>is  non-negative  (positive)  definite  for  all  -»<u«» 


It  is  interesting  to  note  that  positive  real  ness  is  related  to 
passivity.  Otto  Brune  proved  in  the  1930's  that  if  a  filter  (electrical 
network)  was  to  be  realizable  by  passive  circuit  components  (only  resistors. 
Inductors  and  capacitors),  then  its  driving  point  impedance  had  to  be  strictly 
positive  real.  The  resulting  filter  is  a  passive  network,  is  energy  dissipa¬ 
ting,  and  therefore  always  stable. 

2.1.3  Testing  For  Positivity 

In  order  to  apply  the  positivity  theorem  to  establish  the  stability 
of  actively  controlled  LSS,  we  need  a  method  to  test  when  a  square  transfer 
matrix  is  (strictly)  positive  real.  A  time  domain  and  frequency  domain  test 
have  been  found  for  this  purpose. 

Time  Domain  Positivity  Test  T21  Results  obtained  by  B.D.O.  Anderson 
yield  a  relatively  simple  positivity  test  in  the  time  domain.  This  is 
described  next. 

Let  G(s)  be  a  square  matrix  of  rational  transfer  functions  such  that 
G(j»)is  finite  and  G  has  poles  which  lie  in  Re[s]  <  0  or  are  simple  on 
Re[s]  =  0.  Let.  |  A,  B,  C,  G(»)|  be  a  minimal  realization  of  G,  i.e.,  define 


where 


G(s)  =  C(sl-A)”1  B  +  G(-) 


x  =  Ax+Bu 


y  =  Cx+G(®) 

and  A  is  of  least  possible  dimension.  Then  G(s)  is  positive  real  if  and 
only  if  there  exists  a  symmetric  positive  definite  matrix  P  and  matrices 


w  V. 


W  and  L  such  that 
c 


(i)  PA  +  AT?  =  -LLT 
(if)  -6  (®)  +  G^(®) 

(iii)  CT  =  PB  +  LWQ 

The  fact  that  this  test  allows  only  simple  (non-repeated)  poies  on 
the  jco-axis  need  not  worry  us,  since  any  real  structure  has  some  damping,  no 
matter  how  small,  so  that  repeated  modal  frequencies  are  removed  from  the 
jw-axis  by  this  limiting  argument.  Furthermore,  examination  of  Anderson's 
proof  reveals  that  it  would  allow  completely  undamped  modal  frequencies  of 
multiplicity  greater  than  one  because  the  normal  coordinate  formulation  of 
the  structure  allows  strict  diagonal ization  even  with  repeated  eigenvalues. 
The  direct  sum  method  of  modal  subspaces  applies  then  to  each  mode,  regard¬ 
less  whether  the  modal  frequency  is  repeated  or  not. 

G(s)  is  strictly  positive  real  if  for  a  sufficiently  small  real 
o>  0,  G(s)  t  G(s-a)  is  positive  real. 


Note  that  this  statement  is  equivalent  to  modifying  condition  (i) 
of  the  test  to 

PA  +  ATP  =  -Q 

where  Q  =  LlJ  must  be  positive  definite. 


Frequency  Domain  Positivity  Test  [Appendix  of  Reference  1].  A  fre¬ 
quency  test  results  when  one  directly  applies  the  definition  of  positive  real 


To  this  end,  define  the  positivity  index  <5(u.>)  of  a  square  transfer  matrix 
G(s)  by 


<5 (oj )  =  A  m^n  |^[G(jaj)  +  G*(ju))]| 
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where  ^  rnjn  /• l  denotes  "minimum  eigenvalue  of" 

Then 

6 (w )  >  0  for  dll  to  implies  G(s)  is  strictly  positive  real 
<5(w)  -  0  for  all  a)  implies  G(s)  positive  real 
<$( w)  -  0  for  some  w  implies  G{s)  is  not  positive 

A  very  important  advantage  of  this  frequency  domain  positivity  test  over  the 
time  domain  test  is  that  frequency  response  data  may  be  obtained  experimentally 
through  ground  and/or  on-orbit  testing.  This  significantly  reduces  the 
importance  of  detailed  analytical  modelling  of  a  structure  by  finite  element 
methods.  Experimental  frequency  data,  by  necessity,  also  Includes 
actuator  and  sensor  dynamics.  Finally,  it  reduces  the  importance  of  how  to 
model  damping  in  a  structure  as  the  effects  of  damping  are  directly  contained 
in  the  frequency  response  data. 

Positivity  Tests  As  Applied  To  LSS  In  Normal  Coordinates  [1,3] 

Consider  the  use  of  normal  modal  coordinates  in  describing  a  linear  structure. 
The  description  for  G(s)  in  (1)  results  then  in  state  equations  with  the 
system  matrix  ("A-matrix")  in  block  diagonal  form,*  i.e., 


t’h 

where  i  ■+  »  (arbitrarily  large),  u .  is  the  frequency  of  the  j  mode, 

J 

a  *  w.,  and  ?.  is  the  damping  ratio  of  the  jth  mode. 

0  3  0 


Anderson's  positivity  test  can  in  this  case  be  solved  by  hand  for 
arbitarily  large  ;  because  of  the  block-diagonal  form.  An  example  solution 
satisfying  Anderson's  positivity  test  is 


C 


T 


=  B 


*  Assuming  judicious  ordering  of  the  states. 
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This  means  that  the  transfer  matrix  of  a  structure  is  positive  real  if  ideal 
collocated  actuators  and  rate  sensors  are  used.  (But  other  solutions  exist 
also).  Note  that  this  result  is  completely  independent  of  numerical  modal 
data  or  the  number  of  modes  in  the  model.  Note  also  that  application  of  the 
positivity  stability  theorem  to  this  case  demands  that  the  controller 
transfer  matrix  H(s)  be  strictly  positive  real. 


The  inclusion  of  actuator  and  sensor  dynamics  with  the  LSS 
destroys  the  convenient  block  diagonal  form  of  the  modal  state  equations 
and  the  time  domain  positivity  test  (Anderson)  becomes  too  difficult  to 
apply.  The  frequency  domain  test  is. in  this  case,  a  more  attractive 
method  to  test  the  positivity  of  the  infinite-dimensional  plant  (LSS). 

This  is  mainly  because  a  relatively  simple  closed  form  expression  for  the 
transfer  matrix  G0(s)  for  the  structure  can  be  found,  thereby  avoiding  the 
requirement  for  repeated  matrix  inversions  in  the  evaluation  of  the  positi¬ 
vity  index  6(w).  The  terms  appearing  in  the  closed  form  expression  for 
G  (s)  are  given  by 


G(s)  =  Gs(s)  Gq(s)  Ga(s)  and  GQ(s)  = 


where  G(s)  models  the  total  dynamics,  G$(s)  models  the  sensor  dynamics, 
G0 ( s )  models  the  structure,  and  G  (s)  models  the  actuators.  G0(s)  is 

d 

further  defined  by: 


*  +V  +  “j 


■Vv.* 


M*_> 

J  5* 


where  $(r  )  is  the  j  mode  snape  at  the  nr  actuator  location 

ulH 

(z  =  sensor  location), 
s 

Note  that  with  actuator  and  sensor  dynamics  included,  the  positivity 
index  6(<o)  will  now  have  negative  peaks  at  some  modal  frequencies.  This 
indicates  that  the  system  is  no  longer  positive  and  that  the  positivity 
theorem  must  be  augmented  with  operator  embeddina.  This  will  be  the  subject 
?  1  £ 

U  I  v>CV  W  IWtl  t_  *  *  •  V  • 
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There  are  two  basic  steps  involved  when  one  uses  positivity  theory 
for  control  system  design.  One  is  to  assure  that  the  plant  (LSS)  is 
positive.*  The  other  is  to  design  a  controller  that  is  strictly  positive. 
For  the  purposes  of  this  section,  assume  that  ideal  collocated  actuators 
and  rate  sensors  are  used.  This  means,  from  the  arguments  used  in  Section 
2.1.3,  that  the  LSS  is  already  positive.  Hence,  wnat  remains  is  the  design 
of  the  strictly  positive  controller. 


Time  Domain  Controller  Deri<]n  .  Consider  the  linear  quadratic 
Gaussian  (LQG)  controller  topology  shown  in  Figure  2.2.  It  consists  of  a 
state  estimator  followed  by  a  state  feedback  gain.  The  state  equations  for 
the  controller  (without  the  Bc-loop)  are  given  by 


ScCc>  xc  +  V 


If  this  controller  acknowledges  only  the  existence  of  the  controlled 
modes  and  no  special  precautions  are  taken,  the  system  will  surely  suffer 
from  the  effects  of  observation  and  control  spillover.  The  application  or 
the  positivity  design  approach,  however,  will  assure  stability. 

The  basic  approach  is  to  design  the  estimator  gain  matrix  Gc  first, 
and  then  apply  Anderson's  positivity  criterion  so  that  the  transfer  matrix 
H(s)  of  the  controller  is  strictly  positive,  where  H(s)  is  given  by 

*  Although  this  sectiv...  will  presently  assume  the  ideal  situation  where  the 
LSS  is  positive,  Section  2.1.6  will  show  that  the  statement:  "the  plant  is 
positive",  may  be  replaced  by  the  statement:  "the  embedded  plant  is  positive". 
This  generalization  will  permit  the  treatment  of  non-positive  LSS  with  the 
methods  of  this  section. 
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Actuators  Flexlbl 


H(s)  =  Kc  (sl-r)-1  Gc  with  r  -  Ac  -Gc  Cc 

The  estimator  gain  matrix  Gc  is,  for  example,  determined  as  the  stead-- 
state  Kalman  filter  gain  solution,  or  it  is  based  on  pole  assignment,  (i.e. 
a  luenberger  observer),  or  it  is  based  on  any  other  method  that  appears  desira¬ 
ble.  The  feedback  gain  matrix  K.  is  now  observed  to  be  the  output  matrix  of 
the  filter.  So,  from  the  Anderson  test,  if  one  solves  the  Lyapunov  equation 

T  T  A 

p  r  +  r  =  -ll1  •  -  Q 

for  P,  where  Q  is  any  positive  definite  s>mmetric  matrix,  Kc  obtained  from 

K  T  =  P  G 
c  c 

assures  that  H(s)  is  strictly  positive  real. 

Note  here  that  the  choice  of  the  matrix  Q  provides  the  needed  design  freedom 
to  adjust  performance. 


The  Bc  compensation  loop  in  the  state  estimator  is  dashed,  because 
it  cannot  be  explicitly  present  when  K  is  directly  determined  via  positivity. 

V 

This  is  because  the  Bc  loop  mabes  Kc  also  a  part  of  the  system  matrix  of  the 
filter,  rather  than  K,  being  only  an  output  matrix.  But  when  the  state 
estimator  gain  matrix  Sc  -,'s  chosen  to  be  equal  to  the  plant  Inout  matrix  Bc 
associated  with  the  controlled  modes,  I.e.,  G  =  B  -then  the  dashed  loop  may 

U  V 

be  added  after  Kc  has  been  determined  without  violating  che  positivity  of 
the  controller. 


2.1.5  Relationship  Of  The  Time  Domain  Positivity  Design  To  The  General 


With  the  Bc-lcop  included  in  Figure  2-2,  the  topology  of  the  LQG 
controller  a^i  time  domain  positivity  designed  controller  are  the  same.  The 
two  design  approaches  can  result  in  the  same  controller  for  special 
cases,  although  they  are  generally  different. 


In  either  the  LQG  or  positivity  case  the  designer  must  start  by 
selecting  certain  positive  definite  matrices  Q  and  R  and  then,  respectively, 
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solve  Ricatti  and  Liapunov  equations.  If  the  LSS  is  assumed  to  be  described 
in  normal  modal  coordinates  and  collocated  actuators  and  rate  sensors  are 
used.  Table  2.1  shows  a  set  of  conditions  on  the  various  Q's  and  R's  that 
result  in  the  same  controller  for  either  design  approach.* 

The  subscript  "F"  in  Table  2.1  is  associated  with  the  Kalman  filter 
in  the  controller,  and  the  subscript  "OR"  is  associated  with  the  optimal 
regulator  solution  for  K  .  We  define  also  trie  notation 


Qp  3  diag 


with  Q, 


r°  °  1 

.  *  L  .  2  j  =  1.2, ...t 

FJ  L°  4aj°  p  J  C 


=#control led  modes 


where 

varia 


2 

re  'x  s  s.m.  ~  damping  and  c_  any  positive  number  (in  effect  beina 
j  J  J  *  W  ^ 

iance  of  the  observation  noise).  Finally,  Q  is  defined  as 


'0 

0  * 

M 

with  Q.  = 

J 

0 

4a,. 

Tatle  2.1  One  Set  Of  Conditions  For  Which  The  LQG  Solution  Equals 
The  Time  Domain  Positivity  Design. 


K-FILTER/REDUCED  STATE 

DETERMINATION  OF  FEEDBACK 

ESTIMATOR 

GAIN  MATRIX 

LQG 

rf  ■  4 1 

R0R  =  aOR  1 

Qf  .  Qf  .  4  BcBcT 

%R  =  °0R  ^  +  BCBC  * 

POSITIVITY 

SAME  AS  FOR  LQG 

Q  3  Q  ♦  ?BcBcr 

unique,  but  is  one  possible  so’ution 


Relationship  Of  The  Time  Domain  Positivity  Design  To  Other  LSS 
Design  Approaches.  Model  error  sensitivity  suppression  (MESS)  [4]  is  a 
specialized  case  of  LQG  design.  MESS  is  therefore  related  to  the  time  domain 
positivity  design  in  a  manner  similar  to  what  was  just  discussed  previously. 


MESS  with  sensor  throughput  [4]  (shown  in  Figure  2.3a)  is  related 
to  the  positivity  design  since  it  attempts  to  make  the  controller  more 
positive.  That  is. 


H(s)  =  H(s)  +  DI  ,  D  >  C 

+  t  t 

MESS  WITH  MESS  SENSOR 
SENSOR  THROUGH-  THROUGH¬ 
PUT  PUT 

makes  h(s)  more  positive. 

The  direct  "Innovations  feedthrough"  used  by  BALAS  [5]  (shown  in 
Figure  2.3b)  is  a  similar  LQG  type  concept  and  also  results  in  a  more 
positive  controller.  If  either  MESS  with  sensor  throughput  or  LQG  with 
innovations  feedthrough  techniques  are  used  with  collocated  actuators  and 
sensors,  as  various  researchers  have  done,  a  positive  controller  has  in 
effect  been  generated. 


Figure  2.3a  SENSOR  THROUGH-PUT 

Figure  ?.3b  jNNOyATIONS  FEEDTHROUGH 
Figure  2.3  LQG  Variations  Resulting  In  A  More  Positive  Controller 


O  1  ^ 

c  -  I  c 
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Operator  Embedding 


Motivation  For  Using  Embedding.  The  motivation  for  using 
operator  embedding  is  that  strict  positivity  is  too  restrictive  a  require¬ 
ment  to  place  on  the  physical  plant.  Using  Anderson's  results  for  establish¬ 
ing  positivity,  the  transfer  matrix  representing  a  structure  was  shown  to  be 
always  positive  if  ideal  collocatec  actuators  and  rate  sensors  are  used. 
However,  when  the  effects  of  finite  bandwidth  in  the  actuators  and  sensors 
are  included  in  the  dynamic  model  of  the  structure,  its  transfer  function 
ceases  to  be  positive  because  of  the  additional  phase  shift  introduced  by 
these  components.  The  use  o*  pure  position  measurement  and/or  non-col located 
actuators  and  sensors  similarly  destroy  positivity  b>  de-coordinating  the 
spatial  and  time  phase  relationships  existing  between  actuators  and  sensors. 

Strict  positivity  also  limits  controller  design.  This  is  because 
the  difference  between  the  number  of  poles  and  zeros  in  a  positive  transfer 
function  cannot  differ  by  more  than  one.  The  closed  loop  input-output 
control  transmission  roll-off  of  a  design  is  therefore  constrained  to  be  10, 
or  20  dB/decade  for  a  feedback  or  cascade  controller,  respectively.  The 
net  effect  is  that  in  feedback  compensation,  noise  rejection  is  limited, 
while  in  cascade  compensation,  "Bode"  sensitivity  suffers. 

Definition,  Types,  and  Implications  Of  Embedding.  In  order  to 
make  the  positivity  theorem  less  restrictive  and  therefore  more  useful, 
operator  embedding  is  used.  Embedding  essentially  is  a  block  diagram  trans¬ 
formation  of  the  original  system  into  a  stability  equivalent  system*;  one 
then  requires  positivity  of  this  new  system. 

Two  types  of  embedding  transformations  nave  been  considered  in 
this  study.  "F"  embedding,  Figure  below,  refers  to  the  cascade  transforma¬ 
tion,  directly  related  to  the  Popov  multipliers  (simplest  form  is  l+jqu>). 


*  It  is  important  to  note  that  stability  equivalent  or  invariance  does  not 
mean  input-output  transmission  invariance. 


where  F  and  D  are  any  well-behaved  transfer  matrices  (well-behaved  operators) 
including  constant  matrices. 


The  judicious  selection  of  F  and  D  in  the  embedding  process  will 
make  it  possible  to  establish  positivity  of  6  even  though  the  plant  G  is 
not  a  positive  operator.  The  controller  H,  however,  must  "pay"  for  this, 
since  it  must  be  very  stable  and  positive  to  start  with  in  order  that  the 
new  H  will  remain  positive.  To  see  this,  set  F  =  I  and  let  D  be  a  constant, 
diagonal  positive  definite  matrix.  This  clearly  helps  to  improve  the 

positivity  of  G  (think  of  the  trace  of  a  matrix)  but  it  wraps  a  positive 

% 

feedback  loop  around  H,  so  that  H  had  to  be  very  stable  for  M  to  be  (strictly) 
positive  real. 

It  is  important  to  note  that  the  embedding  transformation  is  only 
mathematical  and  is  not  physically  implemented.  The  embedding  may,  however, 
impose  conditions  on  the  controller  which  must  be  physically  implemented.. 

Generally,  embedding  requires  some  knowledge  of  the  plant  to  be 
controlled.  For  example,  constant  "D"  embedding  will  require  that  one  has 
estimated  (albeit,  crudely)  trie  norm  of  the  mode  shapes  and  damping  ratio  of 
the  high  frequency  modes.  If  one  has  additional  knowledge  of  the  plant,  more 
complex  embedding  can  be  used  to  yield  less  conservative  stability  conditions. 

In  a  gross  manner,  embedding  permits  one  to  quantify  stability  conditions  as 
a  function  of  system  knowledge. 

Selection  and  Effects  of  "3"  Embedding.  When  the  plant  G.'s)  is 
nov  positive,  "D"  embedding  may  be  used  to  establish  positivity.  This  is 
done  by  finding  an  operator  D  such  that  3(s)  +  D  is  positive.  The  next  question 
to  be  answered  is  how  D  should  be  chosen  to  assure  that  G"  =  G  +  D  is  positive. 
One  suspects  that  one  must  first  determine  how  negative  G  really  is.  For  this 
purpose  a  positivity/negativity  index  of  G  must  be  determined. 


The  positivity  index  for  G(ju)  was  defined  in  Equation  (2)  as 

=  >  nl-n  {  1  /2[G(  ju)  +  G*  ( jw)]  |  ,  --e[0  ,  ®) 
where  G(j  )  includes  the  structure  plus  actuators  and  sensors. 
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The  operator  0  must  then  be  determined  so  that 

5G+D^  =  ^min{lf  £{G+D)  +  (G+D)*]|>0,  coc[0,®) 

This  last  task  is  considerably  simplified  if  [D(ju)  +  D*  (>)]  is 
chosen  to  have  the  same  eigenvectors  as  [G(jw)  +  G*(ju))].  That  is,  the 
same  eigenvector  matrix  Q  and  its  conjugate  transpose  Q*  (x+x*  is  Hermitian). 

Q*(  jw)[G(joj)  +  G*(jw)]  Q(joj)  ~Ag 
Q*(ja))[D(j(o)  +  D*(j<o)]  Q(jo>)  -  Aq 

where  A^  ^  denotes  "The  matrix  containing  the  eigenvalues  of  (<)“•  For 
this  case,  if  is  chosen  so  that 

<5p(ou)  +  6g  (w)  >0,  V  we  [o,  ®) 

then 

Q*[D  +  D*]  Q  +  Q*[G  +  G*]  Q  >0 
Q*r(o+  G)  +  (D  +  G)*]  Q  >0 

Oq+q  U)  >0  ,  Vwe[o,«  } 

Further  simplifications  exist  if  D(ju>)  =  d  ( jw )  *  I .  For  this  case,  D  (jw) 
always  has  the  correct  eigenvector  structure,  i.e., 

D(j«)  *  d(j«)  •  I  =  d(ju)  .  Q*Q 
so  if  the  real  part  of  d(jy)  is  chosen  as 

Re  [d(jcj]  >  -  »  Vwe  [o,  ®) 

tSor, 
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<Sp+g(o))  >  0,  Vwe  [o,  ®) 


The  effect  of  "D"  embedding  is  to  cut  the  gain  of  the  controller. 

To  see  this,  consider  "D"  embedding  in  Figure  2.1  and  let  H(s)  =  K  >0,  then 
H"  =  (I  -  KD) >0  implies  (I  -  KD)  is  positive. 

Hence 

tr [I  -  KD]  =  tr  [I]  -  tr  [KD]  >0  =  m  >  tr  [KD] 


Where  m  is  the  number  of  inputs  to  the  system.  For  the  case  m  =  1 ,  m  >  tr 
[KD]  =>  k  <  1/D.  Hence  the  gain  is  limited  by  1/D.  This  has  the  obvious 
implications  that  "D"  embedding  gain  stabilizes  the  system.  It  therefore 
places  a  limit  on  the  achievable  damping  and  disturbance  rejection  at  high 
frequencies  that  can  be  achieved  in  a  design*. 


“D"  Embedding  Design  Via  The  Anderson  Test  .  The  previously 
described  time  domain  design  using  the  Anderson  Test  (Page  2-4)  may  be 
used  even  with  a  non-positive  plant.  The  key  difference  when  the  plant  is 
non-positive  is  that  the  minimum  of  the  positivity  index 


<5  .  =  min 
min  w 


mi 


„{?  [G(*.)  +  G*  . 


we 


[0,") 


must  be  computed.  Then  define  a  lSc^"5m-jn  >0  ar,c*  wraP  3  feedback  .■ 

around  the  selected  positive  controller.  This  pre-feedbacK  is  implemented 
and  it  creates  a  new  controller  that  is  Dositive  even  after  the  "D"  embedding 
wraps  a  +6QI  feedback  around  it,  see  Figure  2.4. 


*  It  should  be  noted  that  these  limitations  can  be  mitigated  through 
pre-compensation  of  the  plant  or  through  "F"  embedding.  This  will  be 
shown  later. 
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EMBEDDED  CONTROLLER  H 


Figure  2.4  Embedded  Controller  H  =  H  After  Pre-Feedback  Compensation 

The  design  steps  using  this  approach  are  summarized  as  follows: 

1)  Determine  the  minimum  of  the  positivity  Index,  6ml-n,  of  the  plant 
and  actuators  +  sensors.  Choose  a  >-  <5^ 

2)  Design  a  state  estimator  for  the  truncated  plant,  e.g., 
steady  state  Kalman  filter  or  Luenberger  observer. 

3)  Solve  the  Lyapunov  equation  on  page  2-4  to  obtain  the 
controller  gains.  This  step  will  ensure  that  the  controller, 

H,  Is  positive 

4)  Wrap  negative  feedback  D  *  5QI  around  H  and  implement  this  new 
controller  so  that  it  remains  positive  even  after  embedding. 

The  above  steps  will  ensure  that  the  closed  loop  system  is  stable  by  the 
positivity  theorem,  and  will  result  in  the  controller  topology  shown  in  Figure  2.5. 
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Pre-Feedback  Compensated  Time  Domain  Design 


"D"  Embedding  Design  Using  The  Positivity  (P)-Plot.  "D"  embedding 
may  be  performed  in  a  more  direct  manner  in  the  frequency  domain.  Just  as 
for  the  method  using  the  Anderson  Test,  one  needr  to  first  determine  the 
positivity  index  <s(a>)  of  the  plant  and  then  determine  6Q  >-  6^.  However, 
steps  two  through  four  are  combined  by  predesigning  a  controller  H(s)  to 
withstand  the  "D"  embedding,  i.e.,  such  that 


l 


is 


positive 


To  do  this,  one  uses  a  conformal  map,  called  the  positivity  (P)-Plot.  See 
Figure  2.6.  P-Plot  describes  the  closed-loop  gain  and  phase  of  a  positive 
feedback  system  as  a  function  of  its  open-loop  values.  Hence,  all  one  has 
to  do  is  make  sure  that  the  controller  K(s)  with  loop  gain  &  ,  lies  within 
the  -  90°  closed-loop  phase  region  of  the  P-Plot  (shaded  region). 

"D"  Embedding  Design  Combined  With  Frequency  Domain  Prccomponsa- 
tion.  In  practice,  precompensation  of  the  plant  may  be  needed  in  order  to 
meet  stability  and  performance  objectives.  One  can  combine  "M"  precompensa¬ 
tion  and  "D"  embedding  by  first  carrying  out  a  full  multivariable  frequency 
domain  type  design  (summarized  next).  This  results  in  the*  determination  of 
a  cascade  precompensator  M($).  The  net  effect  is  a  cascade  pr  .*-ca.ipen sated 
system  within  which  the  controller,  M(s),  has  already  bjen  designed  for  per¬ 
formance  (low  sensitivity,  noise  rejection  and  damping)  and  desired  high 
frequency  roll-off.  The  latter  is  important  since  the  maximization  of  the 
roll-off  minimizes  uncertainties  in  the  "D"  embedding  required.  Once  the 
performance  design  is  completed,  one  can  determine  the  "D"  embedding  operator 
that  makes  the  precompensated  plant  positive.  This  is  typically  done  by 
determining  the  positivity  index  s(u)  of  the  precompensated  plant,  finding 
its  most  negative  peak  6  ^n,  and  choosing  D  =  I ,  SQ  >-  >0.  The  "H" 

part  of  the  controller  consists  then  merely  of  gains  which  must  be  set  to 
less  than  or  equal  to  1 /6q .  The  actually  Implemented  controller  satisfying 
the  positivity  design  is,  as  shown  in  Figure  2.7,  then  given  by  (1/6  )M(s). 
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PHRSE  (OF.GRFJ  S  ) 

FIGUr.E  2.6  Positivity  (P)  -  Plot 


CONTROLLER 


PRE- COMPENSATED 
PLANT 


Figure  2.7.  "D"  Embedding  A  Precompensated  Plant. 


The  key  to  the  extension  of  the  classical  Nyquist/Bode  approach  to 
the  multivariable  case  [6]  are  the  characteristic  gains  of  the  transfer  matrix, 
G(s),  describing  the  open-loop  plant.  The  characteristic  gains  A  (s)  are 
defined  as  the  solutions  to  the  equation 

GET  [A  (s)  •  I  -  fi(s)]  =  0  ;  se  [oj,  -j) 

The  importance  of  the  characteristic  gains  stem  from  the  fact  that  the 
characteristic  gains  of  a  unity  feedback  system  are  related  to  the  open-loop 
values  in  a  manner  analogous  to  single-input-single-output  systems,  i.e.. 


closed-loop  gain 


open-loop  gain 
1  +  open-loop  ga in 


The  main  results  using  this  technique  are: 

1)  Closed-loop  stability  occurs  if  and  only  if  the  net  sum  o * 
counter-clockwise  encirclements  of  the  (-1,0.)  point  by  the 
characteristic  gains  is  equal  to  the  number  of  open-loop 
unstable  poles. 

2)  To  achieve  low  control  interaction  in  a  frequency  band,  it 
is  sufficient  that  all  the  character!' stiu  gains  have  a  large 
magnitude  over  that  frequency  band. 
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3)  In  frequency  bands  of  'low  interaction,  the  notions  of  gain 
margin  and  phase  margin  may  be  applied  as  a  qualitative 
assessment  of  performance. 

4}  I  here  is  a  relationship  between  tracking  nCCui'aCy  aiid  the 

magnitude  of  the  characteristic  gains.  Under  the  appropriate 
assumptions  the  relationship  is  analogous  to  the  single-input- 
single-output  case. 


2.1.7  General  Formats  For  Computer  Aided  Design 

The  functions  to  be  evaluated  in  the  positivity  design  are  lengthy 
and  virtually  impossible  to  do  by  hand.  A  computer  aided  design  approach 
is  therefore  required.  Figures  2.8a  and  2.8b  respectively  show  the  general 
software  formats  for  computer-aided  time  domain  design  and  frequency  domain 
design.  Enclosed  by  the  parenthesis  are  listed  TRW  programs/sub-programs 
which  implement  and  perform  the  required  computations.  The  key  elements 
required  to  make  the  package  useful  are  a  fast  computer,  highly  accurate 
linear  al  oebra/eigerivalue  subroutines,  and  a  highly  interactive  man-machine 
iterative  loop  with  graphics.  The  current  TRW  package  employes  this 
approach  and  makes  extensive  use  of  state-of-the-art  subprograms  and  graphics 
from  the  IMSLIB,  EISPAK  and  TGSLIB  Libraries. 


I - 1 

-H  ! 


-  4  ON-LINE  i 
HUMAN  t 
J  INTERACTION  | 


Figure  2.8a  Computer  Programs  Required  For  Time-Domain  Design 
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FLAT  PLATE  E*  ?LES 


As  an  example  consider  a  simply  supported  (pivoted,  4  sides), 

1.5  x  1.5  m  aluminum  plate  that  is  1.5mm  thick,  shown  in  Figure  2.9.  The 
first  10  modal  frequencies*  and  a  closed  form  expression  for  the  mode 
shapes  are  given  in  the  Table  2.2.  Assume  the  objective  is  to  actively 
control  the  first  three  modes.  Since  repeated  modal  frequencies  are  present, 
one  requires  at  least  two  sensor/actuator  pairs  to  ensure  observability/ 
controllability.  Two  point-force  actuators  and  two  rate  sensors  collocated 
have  therefore  been  chosen  as  shown  in  the  Figure. 

2.2.1  Idealized  Time  Domain  Design  Example 

In  the  design  assume  we  have  idea1  collocated  actuators  and  r«v 
sensors.  The  plate  transfer  matrix  will  therefore  be  positive  real  regardless 
of  the  exact  values  of  the  modal  frequencies  and  modal  truncation.  The  basic 
time-domain  approach  Is  then  to  design  a  strictly  positive  real  controller 
transfer  matrix  H(s)  by  applying  Anderson's  positivity  criterion. 

Referring  to  Figure  2.2,  first,  since  three  modes  are  to  be  con¬ 
trolled,  a  sixth  order  state  estimator  Is  designed  by  solving  for  the  steady- 
state  Kalman  filter  gain,  Gc,  of  the  truncated  model  (3  controlled  modes). 

Then  one  notes  that  K  is  the  output  matrix  of  the  controller,  i.e., 

V 

H(s)  =  K  (Sl-r)"1  G  with  rs  A  -G  C 
c  c  c  c  c 

So,  from  Anderson's  test,  one  can  solve  a  Lyapunov  equation  for  the  state 
feedback  gain  Kc  which  will  ensure  that  the  controller  is  positive.  LLT=Q= 
Diagonal  (100),  here,  was  found  to  yield  good  performance.  Since  the  plate 
Is  positive  and  the  controller  is  designed  to  be  strictly  positive,  the 
design  is  guaranteed  to  be  stable  even  when  an  arbitrary  number  of  residual 
modes  are  retained. 


*  Damping  is  .4% 
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Table  2-2  Modal  Parameters  For  Plate  Example 


Wo4a  Sfeapa  Wilt! 

Ikaktr  Intagars  Fraquatic}' 

_ i _ kl _ Hx 


1 

l 

1 

3.24 

2 

i 

2 

8.10 

3 

2 

1 

8.10 

4 

2 

2 

12.  H 

5 

3 

1 

18.20 

6 

1 

3 

16.20 

7 

2 

3 

21.00 

• 

3 

2 

21.00 

9 

4 

1 

27.3 

10 

1 

K 

27.3 

Dm  m4i  ikif 

m  nr* 

flttlllli 

Iran  tha  aquation 

♦jfc't  y)  • 

.449  ala 

(  k,t 

1  1.3 

-) 

Figure  2.9  Flat  Plate  With  Collocated  Actuators/Senscrs 


2-27 


Figures  2.10a-d  show  closed-loop  simulation  results  for  a  plate  model 
retaining  the  first  10  modes.  Figure  2.10a  and  b  show  the  plate  deflection 
transient  response  and  control  effort  at  the  control  stations.  Figures  2.10c 
and  d  show  that  the  transient  response  of  the  controlled  modes  Is  good  and 
the  residual  modes  are  stable. 

2.2.2  Time  Domain  Example  With  ::D"  Embedding 

Using  the  same  flat  plate  model,  the  "D"  embedding  design  method 
using  the  Anderson  test  will  be  used  to  demonstrate  the  positivity  time- 
domain  technique  with  finite  bandwidth  actuators  and  sensors.  The  actuators 
and  sensors  have  been  modelled  as  first  order  lags  with  a  bandwidth  of  30  Hz. 

Following  the  design  steps  on  page  2-17  one  first  computes  the 
positivity  index  <Sg(oi).  This  is  plotted  in  Figure  2.11.  Here,  one  notes 
that  because  of  the  collocated  actuators  and  sensors  6^)  is  positive 
except  for  some  negativity  at  the  modal  frequencies.  This  effect  is  due  to 
the  phase  shift  in  the  actuators  and  sensors. 

Next,  as  in  the  first  example,  one  computes  the  Kalman  gains  6 

c 

for  the  3-control led-mode/ truncated  plate  mode  and  computes  the  state 
feedback  gain  K  according  to  the  Lyapunov  equation.  This  guarantees  that 

V 

H(s)  is  positive. 

Finally,  one  compensates  H^s)  so  that  one  can  embed  the 
system  and  have  Sq+qU)  >0  for  all  u.  The  simplest  way  to  do  this  Is,  as 
described  earlier,  to  choose  D  equal  to  5q  -I,  where  6q  >  >0  and  6m1n 

is  the  most  negative  value  of  6q(u).  =  0.15  is  a  good  choice  here. 

The  implemented  controller 

Hj  •*  (I  +  HD)-1  H 

will  therefore  permit  a  fictitious  "D"  embedding  with  D  =  0.15.1  and  still  remain 
positive.  This  leaves  the  closed-loop  system  looking  like  that  shown  in 
Figure  2.5  and  it  is  stable  by  positivity  theory. 
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Figures  2.12a  and  b  show  the  simulated  transient  response  of  the 
designed  system.  Figure  2.12a  shows  the  transient  Is  stable.  Figure  2.12b 
shows  that  the  control  effort  Is  very  modest.  It  Is  Interesting  to  compare 
these  simulation  results  with  those  of  the  Idealized  system,  figure  2.10a-d. 
The  transient  In  the  present  design  Is  half  as  fast  and  the  control  effort 
Is  reduced  by  a  factor  of  50.  These  results  make  good  Intuitive  sense  since 
one  Is  now  dealing  with  finite  bandwidth  actuators  and  sensors. 


2.2.3  Frequency  Domain  Precompensation  Example 

Using  the  same  flat  plat  as  an  example,  the  combined  precompensatlcns 
and  "D"  embedding  technique  outlined  on  page  2-20  will  be  used  to  show 
how  one  can  manipulate  the  characteristic  gains  and  the  positivity  Index. 

The  characteristic  gains  for  the  plate,  shown  In  Figure  2.13a,  lead 
one  to  believe  that  the  system  has  sufficient  stability  margins  to  permit 
high  loop  gains.  However,  the  positivity  Index,  shown  In  Figure  2.13b  has  not 
decreased  at  high  frequencies,  a  region  where  high  uncertainty  exists  in 
the  modal  data.  Hence,  one  needs  to  use  low  loop  gains  In  order  to  accommo¬ 
date  for  a  conservative  estimate  of  the  SQ  used  during  the  "D"  embedding 
operation.  This  of  course  is  at  the  expense  of  performance.  To  avoid  this 
problem,  one  would  like  to  shape  the  loci  so  that  a  decreasing  positivity 
index  is  obtained  at  high  frequencies.  If  this  is  done  such  that  there  is 
little  or  no  effect  on  the  controlled  modes,  it  would  permit  a  less  conserva¬ 
tive  estimate  of  <SQ,  the  latter  permitting  higher  loop  gains.  The  approach 

that  has  been  used  is  to  shape  the  characteristic  gains  and  the  positivity 
index  by  using  a  commutative  compensator  [6], 
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The  Figures  2.14a  and  b  show  the  effect  of  4th  order  compensation. 

This  compensation  was  optimized  for  attenuation  of  the  last  two  modes  and 
the  final  "cross-over",  and  for  minimum  phase  lag  at  the  controlled  mode 
frequencies.  The  compensai.  '  characteristic  gains  and  positivity  index 
show  that  this  objective  was  met.  The  index  at  the  highest  frequency  has 
been  reduced  by  a  factor  of  three.  Thus  permitting  a  less  conservative 
estimate  of  to  be  used  without  risking  instability  at  some  unmodeled  high 
frequency  mode. 

2.3  DRAPER  MODEL  #2  EXAMPLE 

In  this  section,  the  Draper  Model  #2*,  shown  in  Figure  2.15a  will 
be  used  as  a  test  evaluation  structure  for  the  positivity  design  approach. 

The  objectives  are  two-fold.  One  is  to  demonstrate  how  the  technique  can 
be  applied  to  establish  the  stability  of  the  structure.  Inclusive  of 
actuator  and  sensor  dynamics  and  rigid  body  control.  The  second  is  to  show 
how  one  can  couple  the  technique  with  on-orbit  tuning  to  create  a  stability 
ensuring  design  methodology  that  can  meet  stringent  performance  goals. 

The  np,-formance  goals  for  Draper  Model  #2  are  to  maintain  the 
1 ine-of-sight  (LOS)  rotation  and  DEFOCUS  errors  below  .05  u  rad  and  25  u  m, 
respectively.  This  accuracy  is  to  be  maintained  even  though  sinusoidal 
disturbances  are  acting  directly  on  the  structure.  Figure  2.15b  shows  the 
locations  and  amplitudes  of  these  disturbances. 

Conceptual  Design.  One  way  to  understand  the  LSS  control  problem 
is  to  think  of  the  LSS  and  the  surrounding  environment  as  a  closed  system 
wherein  energy  is  stored,  transformed  and  dissipated.  The  vibration 
suppression  problem  is,  in  these  terms,  to  simultaneously  maximize  energy 
dissipation  within  the  "clean"  structure,  and  minimize  the  vibrational 
energy  created  in,  and  transmitted  to  it  from  the  outside.  Ideally,  the  latter 
would  need,  to  he  totally  performance  robust,  to  spatially  isolate  the  whole 
clean  structure  from  the  outside  over  all  energy  frequencies.  This  cannot  be  done  in 

*  A  description  of  Draper  Model  #2  is  given  in  Appendix  A 
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2.15a  DRAPER  MODEL  #2 
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'•.eneral,  however,  because  total  active  spatial  isolation  would  need  to  be 
continuously  distributed  and  have  infinite  bandwidth.  Fortunately,  in 
practice,  the  most  energetic  sources  are  spatially  lo-alized  and  a  large 
percentage  of  the  emitted  power  is  over  a  soec^fic  band  of  frequencies. 

Hence,  a  practical  solution  to  the  isolation  problem  exists  when  energy 
sources  can  be  identified  as  to  content  and  locale.  A  tradecff  between 
the  amount  of  isolation  of  the  "clean"  structure  from  the  equipment  (dirty) 
section  and  amount  of  dissipation  also  exists,  the  amount  of  dissipation 
required  depending  on  the  energy  "leaked"  by  the  isolating  mechanism. 

The  strategy  for  Oraper  Model  12  will  be  to  put  the  burden  of  the 
prcbles  upon  the  isolators.  This  will  take  advantage  of  the  fact  that 
the  locale  and  energy  content  of  the  disturbance  sources  are  well  known. 

A  passive  Isolator  system  will  therefore  first  be  used  to  isolate  the  clean 
structure  from  the  disturbances  in  the  dirty  box;  then,  spatial  locations 
for  passive  damping  will  be  chosen  so  that  leaked  energy  in  the  disturbance  band 
will  be  maximally  dissipated;  the  rigid  body  controller  will  then  be  designed 
to  provide  attitude  control  and  isolation  to  the  entire  structure  from  low 
frequency  disturbances  emanating  from  the  environment;  finally,  the  clean 
structure  disturbance  will  be  locally  dissipated  and  isolated  from  the 
"spatial  entity"  LOS.  It  is  important  to  note  that  vibration  in  the 
entire  clean  structure  need  not  be  suppressed.  It  is  sufficient  to  suppress 
only  the  spatial  entity  described  by  1 ine-of-sight. 

Physical  Considerations.  As  in  any  other  feedback  control  system, 
the  quality  of  the  actuators  and  sensors  must  be  commensurate  with  the 
desired  closed-loop  performance.  This  leads  to  some  prerequisite;  -hat 
must  be  met  by  the  actuators  and  sensors  if  the  overall  performance  of 
.05  ^rad  error  is  to  be  achievable. 

With  respect  to  the  sensors,  the  desire  to  point  with  LOS  errors 
less  than  .05  urad  means  that  one  must  have  at  least  one  inertial  sensor 
capable  of  providing  a  reference  down  to  this  level  Second,  the  proper 
control  of  structural  vibrations  will  require  the  active  control  sensors 
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to  he  sensitive  enough  to  measure  the  expected  deflections  and  rotations 
incurred  at  the  sensing  station  during  a  .05  u  rad  LOS  deviation.  Finally, 
the  disturbance  rejection  requirenents  imply  that  the  active  control  sensor 
bandwidth  must  exceed  the  disturbance  bandwidth  by  a  factor  of  5-20.  The 
latter  depending  upon  the  type  of  expected  disturbance,  the  bandwidth  of  the 
actuators,  and  the*  mode  shapes  and  spacings  of  the  modal  frequencies  between 
the  control  bandwidth  and  the  actuator  bandwidth.  Similar  bandwidth 
considerations  must  be  given  to  the  actuator  mechanisms.  However,  the 
actuator  must  first  be  capable  of  acting  upon  an  LOS  deviation  of  .05  y  rad. 
This  means  that  its  deadband  must  be  narrower  than  the  expected  deflections 
and  rotations  Incurred  at  the  control  stations  during  the  LOS  deviation. 

The  decision  of  whether  or  not  an  actuator-sensor  pair  should  work 
on  a  relative  or  inertial  basis  also  has  to  be  considered  carefully. 
Intuitively,  an  inertial  disturbance  {ex.,  moving  masses  in  a  cryo-cooler) 
acting  on  a  totally  rigid  object  can  only  be  counteracted  using  inertial 
controls.  Relative  control,  on  the  other  hand,  for  use  in  vibration 
suppression  is  as  likely  to  redirect  energy  as  much  as  dissipate  it.  After 
all,  this  controller  has  its  inputs  connected  to  a  moving  system.  Hence 
the  controller  can  only  dissipate  energy  as  fast  as  its  reference  permits 
it.  Inertial  control  Is,  however,  referenced  to  something  that  is  theoreti¬ 
cally  not  moving,  an  inertial  frame  where  ideally  no  dissipation  rate  limits 
will  exist.  It  seems  reasonable  to  expect,  therefore,  that  inertial  distur¬ 
bances  will  require  a  higher  number  of  relative  controllers  than  inertial 
controllers  for  an  equivalent  rate  of  dissipation. 

Determining  Modes,  Nodes  and  Truss  Members  To  Be  Controlled. 

Previous  esperience  by  TRW  on  the  control  of  flesible  spacecraft  and 
missiles  (ex.,  TORS  and  Minuteman)  has  shown  that  the  G-factor  gives  a 
good  measure  of  the  relative  amplitude  of  a  mode  at  a  particular  node. 

The  G-factor  is  defined  as  the  d.c.  gain  of  the  transfer  function  between 
two  nodes  at  a  modal  frequency  of  interest,  i.e.. 
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i,  j  are  node  members 
n  is  the  mode  number 

k  identifies  the  spatial  coordinate  (e.g.  translation: 

x.  y,  z;  rotation  X  ,  x  ,  x  ;  or  along  a  truss  direction). 

x  y  z 

The  TRW  computer  program  MDLSFT  automates  the  sifting  process.  It 
will  go  through  all  the  modal  data  and  will 


1)  Rank  modes  by  5-factor  given  Important  nodes  or  mode  shapes. 

2)  Rank  nodes  by  G-factor  given  Important  modes. 

3)  Rank  relative  node  members  given  important  modes 

(those  nodes  connected  by  a  truss  member)  by  G-factor  along 
the  line  connecting  the  nodes. 


2.3.1  Control  Topology  For  Rigid  Body  and  Active  Structure  Control 

A  functional  block  diagram  of  the  control  topology  used  on  Draper 
Model  #2  is  shown  in  Figure  2.16,  The  controller  H(s)  is  partitioned  as 
shown  in  Figure  2.17.  Hp(s)  controls  the  flexible  dynamics  and  ( s ) 

controls  the  rigid  body  attitude.  The  idea  here  is  to  make  use  of  the 
fact  that  if  both  H^(s)  and  Hp(s)  are  positive,  then  the  composite 
controller 


Ha<5> 

0 

0 

Ms) 

r 

is  also  positive.  These  controllers  can  therefore  be  designed  separately 
and  because  each  one  is  positive,  either  Hp  of  Hj,  nay  be  set  to  zero  without 
destabilizing  the  system. 
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In  order  to  make  the  plant  as  positive  as  possible,  over  the  largest 
frequency  range,  the  rigid  body  controller  assumes  only  rate  is  measured. 
Attitude  control,  however,  will  require  some  position  feedback.  HR(s)  is 
therefore  chosen  as 

km 

Vs>  ■  kr,  +  — 

Now  HR(s)  is  positive  real  for  positive  semi-definite  and  the 

positivity  theorem  with  embedding  can  be  easily  applied.  Note  here  that 
for  spacecraft  using  integrated  gyro  rate  for  attitude  reference,  H^(s) 
models  the  true  physical  situation. 

2.3.2  The  Positivity  Design 

Following  the  conceptual  design  procedure,  the  approach  consisted 
of:  1)  the  passive  isolator  design;  2)  the  passive  member  damper  design; 

3)  the  rigid  body  design;  and  finally  4)  the  active  damping  design. 

The  first  quantities  that  were  computed  are  the  LOS  disturbance 
rejection  requirements.  This  was  computed  by  normalizing  the  dB  scale 
according  to  OdB**  s  1  rad/N.  Here, .05  p  rad/N  represents  -146  dB,  so  for 
the  2 ON  and  ION  disturbances  considered  in  this  study,  it  can  be  shown  that 
the  requirements  are  equivalent  to: 


@  31 .42  rad/sec 

@  62.83  rad/sec 

From  Top  (Node  3.  ) 

- 

-166  dB 

From  Bottom  (Node  46) 

-172  dB 

- 

Passive  Isolator  Design.  The  passive  isolator  design  is  an  example 
of  the  ideal  positivity  design.  Here,  rate  gain  is  the  implicit  feedback 
parameter,  and,  since  there  are  no  actuator  and  sensor  dynamics,  the  dash- 
pots  create  a  natural  positive  definite  feedback  gain  matrix. 


★ 


Gain  in  dB  will  refer  to  20  LCG,q  [magnitude] 
stated  otherwise. 
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unless  specifically 


V- 


The  approach  used  to  determine  the  actual  Isolator  dashpot  constants 
is  based  on  the  multivariable  frequency  domain  approach  of  Mac  Farlane  [6]. 
For  the  problem  here,  this  requires  that  one  determine  the  characteristic 
gains  and  transmission  zeros*  of  the  system.  The  characteristic  gains  were 
used  as  a  guide  in  the  selection  of  loop  gains  resulting  in  desirable 
sensitivity  and  damping  properties.  Th'  transmission  zeros  were  used  to 
determine  which  eigenvalues  will  become  real. 

If  one  looks  at  the  frequency  response  of  a  normalized  2nd  order 
transfer  function  with  ?=1,  one  will  notice 

PHASE  @  2m  -  PHASE  G>  \  w  =75° 
n  2  n 

Hence,  a  plot  of  the  characteristics  gains  in  the  Nichols  chart  allows  one 

to  approximately  determine  the  gains  g^  for  which  the  eigenvalues  will  first 

touch  the  real  axis.  That  is,  when  the  closed-loop  phase  contours  changes 

between  +38°  and  -38°  as  m  changes  between  Iw  and  2m  . 

c  n  n 

An  eigenvalue  analysis  for  loop  gains  j  g.,  g.  and  2g.  then  permits 
a  quick  check  to  see  if  isolator  damping  is  adequate,  and  also  supplies  a 
selection  of  possible  alternatives.  Final  selection  of  gains  was  based  on 
the  set  providing  the  best  combination  of  overall  damping  and  isolation. 

The  frequency  response  plots  in  Figures  2.18a  and  b  show  the  results, 
of  the  design.  The  responses  are  from  the  disturbance  source  locations  to 
the  1 ine-of-sight  outputs.  The  defocus  error  is  not  shown  as  it  easily 
meets  the  ?5  ^  m  design  goal. 

Passive  Member  Damper  Design.  The  passive  member  damper  design  is 
another  example  of  the  ideal  positivity  design.  Damper  values  were  therefore 
determined  by  using  the  same  procedure  outlined  in  the  passive  isolator 
design.  The  critical  part  was  determining  which  and  how  many  truss  members 


*  Eigenvalues  of  the  system  as  the  feedback  gains  approach  infinity. 
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4)  Due  to  the  large  dimensions  of  the  models  involved,  the  frequency  response  plots  are  only  a 
first  order  approximation  to  the  true  response.  The  approxinw  tion  is  based  on  ignoring  the 
modal  coupling  introduced  by  the  insertion  of  the  isolator. 


were  to  have  the  dashpots  in-line  with  themselves.  In  order  to  determine 
this,  the  MDLSFT  was  run  using  two  options:  1.  Given  Important  mode 
shapes  (LOS),  rank  modes  by  G-factor;  these  resulting  modes  were  then  used 
as  the  Important  modes  in  the  option,  2.  Given  Important  modes,  rank 
relative  node  members  by  G-factor.  The  member  dampers  were  then  chosen 
from  this  list  according  to  the  magnitude  of  their  modal  influence  and  the 
number  of  modes  they  Influenced.  The  number  of  members  was  determined  as 
the  minimum  number  to  have  Influence  on  all  the  critical  inodes. 

Four  member  dampers  resulted  using  the  above  procedure.  The  fre¬ 
quency  response  plots  In  Figure  2.19a  and  b  show  the  results.  These 
figures  Indicate  improvement  In  the  disturbance  rejection  response,  but 
active  control  is  still  required  to  meet  specifications. 

Rigid  Body  (Attitude)  Control  Design.  The  rigid  body  controller 
Is  the  first  place  In  the  Model  #2  design  wherein  active  control  will  be 
required.  The  ACS  is  located  on  node  44  of  the  dirty  box  and  it  is  composed 
of  CMG's,  rate  gyros  and  a  three  axis  controller.  The  CMG's  were  assumed 
to  have  a  5  Hz  bandwidth,  the  gyros  10  Hz.  This  controller  was  designed 
using  the  "D"  embedding  design  procedure  outlined  earlier.  The  main  design 
steps  were  to: 

1)  Compute  the  positivity  index 

6q(“)  =  *  min  {1  +  v  uz  C°»“) 

where 

G(jcu)  is  the  input-output  transfer  matrix  of  the  rigid 
body  of  the  structure  including  actuator  and  sensor  dynamics. 

2)  Determine  a  fictitious  diagonal  operator  D(jo>)  =  I-d(ju) 
such  that  5q(oj )  +  Real  [d(ju>)]  >  0  for  all  u. 

3)  Physically  implement  a  positive  controller  H(ju)  that  remains 
positive  even  after  fictitious  positive  feedback  D(jto)  has  been 
wrapped  around  it. 
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In  the  present  design  the  rigid  body  controller  was  assumed  to  be 


of  the  form 


=  KRA  + 


where  KRA  =  and  KpA  =  I-kDA«  I£  can  be  shown  that  with  tcR)  *  1.2  x 

106  N-m  sec/rad  and  kpA  =  8  x  10^.  N_m/rad,  Hp(s)  will  remain  positive  with 
fictitious  "D"  embedding  tracking  the  -<5„U)  curve.  Figures  2.20a  and  b 

U 

show  the  positivity  index  6q(u)  and  the  Nichols  plot  of  the  characteristic 
gains  for  the  rigid  body.  Also  shown  on  Figure  2.20b  is  the  resulting  rate  loop 
operating  point. 

Active  Dc-aping  System  (Positivity  Design*)-*  The  design  methodology 
for  the  active  damping  system  is  the  same  as  in  the  rigid  body  design.  The 
objective  here,  however,  is  to  maximally  isolate/dissipate  the  energy 
emanating  from  the  disturbance  sources. 

The  active  control  locations  were  chosen  to  be  those  nodes  closest 
to  the  disturbances  (46  and  37).  It  can  be  argued  that  these  locations  may 
be  sensitive  to  disturbance  location,  but,  it  can  also  be  argued  that  they 
are  the  most  robust  w.r.t.  to  the  modal  data.  Placement  In  real  situations 
would  tnerefore  depend  on  trading  the  degree  of  expected  spatial  uncertainty 
for  the  degree  of  expected  modal  data  uncertainty. 

Ideally  one  model  should  be  used  for  determining  the  positivity 
index  of  the  active  control  system.  However,  the  large  number  of  modes 
within  the  actuator/sensor  bandwidth*made  this  difficult.  Instead,  two 
models  were  determined  for  computing  a  worst  case  index.  One  model  was 
based  on  G-factor  rankings  of  those  modes  affecting  node  37  the  most,  and 
the  other  on  those  affecting  node  46  the  most.  The  resulting  positivity 
plot.  Figure  2.21a  shows  the  node  46  model,  which  is  the  most  "negative". 

Rate  gams  and  position  compensation  were  therefore  based  on  this  model. 

Figure  2.21b  shows  the  rate  loop  operating  point  in  the  Nichols 
plot  for  the  active  system  model.  The  position  compensator  was  chosen  as 

*Actuator/sensor  bandwidth  is  100  Hz. 
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th.i  o.-i  i  2nd  order  filter  which  would  simultaneous1/  boost  the  gain  at 
30  and  60  rad/sec  to  +9  dB  and  minimize  high  frequency  gain.  This  choice 
permits  frequency  selective  high  loop  gains  at  the  frequencies  requiring 
the  most  disturbance  rejection.  Higher  order  filters  could  have  been 
chosen  here  for  greater  selectivity,  but  a  better  characterization  of  the 
uncertainty  in  the  disturbances  would  have  been  required.  If  the  distur¬ 
bances  were  truly  sinusoidal  and  at  exactly  the  given  frequencies,  the 
optimal  compensator  would  have  had  one  undamped  pole  for  each  sinusoidal 
disturbance  frequency. 

2.3.3  Results  For  The  Positivity  Design 

Figure  .2.22  shows  a  block  diagram  of  the  entire  Draper  Model  #2 
positivity  designed  system.  This  system  was  analyzed  and  simulated  to 
check  the  robustness  and  performance  of  the  design. 

To  analyze  the  positivity  design  performance,  a  nominal  evaluation 
model  retaining  the  largest  G-factors  in  the  LOS  was  constructed.  The 
frequency  response  of  this  model,  shown  in  Figures  2.23a,  b,  c  was  then 
computed  from  various  points  of  interest.  The  nominal  design  yielded  a 
rigid  body  bandwidth  of  at  least  .4  rad/sec,  and  disturbance  rejections 
from  disturbance  source  to  LOS  as  shown  below.  Clearly,  the  specification 
could  not  quite  be  attained  with  the  untuned  positivity  design.  This  was 
more  or  less  expected,  since  the  positivity  design  is  conservative. 


LOS  x 

LOS 

y 

req'd  to 
meet  spec 

From  top  (node  37)  at  62.83  rad/sec 

-152  dB 

-145 

dB 

-166  dB 

From  bottom  (node  46)  at  31.42  rad/sec 

-170  dB 

-158 

dB 

-172  dB 

The  root-sum-squared  prediction  based  on  this  data  is  LOS  x  error  =  .258  yrad 
and  LOS  y  error  =  .751  yrad.  These  predictions  agree  with  those  obtained 
through  simulation,  as  shown  in  Figures  2.24a,  b,  c.  Figures  2.24d  and  e 
show  that  the  required  actuator  force  is  a  very  reasonable  25N,  which  can  be 

implemented  with  a  linear  momentum  exchange  device(proof  mass  actuator)  weighing 
about  5  Kg.  z~49 
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Figure  2.22  Block  Diagram  of  DRAPER  Model  #2  Design 


Figure  2.23a 

Rigid  Body  Command  To 
Rigid  Body  Sensor 
Frequency  Response. 


Figure  2.23b 

Disturbance  Rejection 
From  Node  37  to  LOS 


Figure  2.23c 

Disturbance  Rejection 
From  Node  46  to  LOS 


NOTES:  Model  retains 
modes  4-7,  21 , 
28,  30,  34,  35, 
37,  39,  40,  41 


i 
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Figures  2.24f-h  show  the  ACS  step-input  transient*Ax ,  A  ,  and  Az, 
were  simultaneously  commanded  .01,  -.01  and  .01  rad,  respectively.  These 
Figures  show  a  stable,  smooth  rise  in  less  than  5  seconds  and  with  little 
ringing  from  the  multitude  of  bending  modes. 

In  order  to  demonstrate  the  robustness  of  the  positivity  design,  the 
nominal  evaluation  model  was  perturbed  as  described  in  Appendix  A.  The 
simulation  results  using  perturbed  Model  #2  and  #4  show  that  the  system  is 
stable.  Moreover,  the  performance  change  from  the  nominal  case  is  imper- 
ceptable*.  See  Figures  2.25a-f. 


2.3.4  Positivity  Design  "Tune-up*1 


The  results  of  the  previous  design  indicate  that  the  disturbance 
tracking  has  to  be  improved  in  order  to  meet  all  the  performance  require¬ 
ments.  Since  the  open  and  closed-loop  characteristic  gains  of  this  design 
are  related  analogously  to  the  Bode/Nyquist  single-input-single-output 
case,  the  improvement  in  tracking  is  realized  by  increasing  the  magnitude  of 
the  characteristic  gains. 


i .e. ,  for  large  a(s) , 


tracking  error 


A  ( s ) 

i  +  (Kp  +  krsTTTs7 


i 

+  V 


S  in  the 
disturbance 
band 


It  was  now  assumed  that  the  characteristic  ains  in  Figures  2.20b 
and  2.21b  resulted  from  ground  processing  of  post  un-orbit  test  data  and 
describes  the  system  accurately  (note  that  the  transfer  matrix  can  be 
obtained  from  on-orbit  frequency  response  data  via  the  FFT).  Therefore 
much  more  elaborate  embedding  can  be  used,  or  if  the  system  can  be  assumed 
to  be  linear,  multivariable  frequency  domain  techniques  can  be  applied. 

The  locus  in  Figure  2.21b  predicts,  then,  that  the  DC  gain  increase  required 
to  meet  specification  would  result  in  an  unstable  system,  so,  either  addi¬ 
tional  compensation  or  increased  actuation/sensing  bandwidth  is  required. 

In  the  interest  of  time,  we  doubled  the  active  system  actuator/sensor 
BW  to  1200  rad/sec  so  that  we  could  assure  that  the  required  1oop  gains 
would  not  destabilize  the  system. 

*  As  measured  by  sterr'y  state  peak-to-peak  error. 
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Figures  2.26  a-c  show  the  LOS  errors  for  the  "tuned-up"  system 
as  simulated  on  the  computer.  The  tune-up  design  meets  all  specifications. 


NOTES:  1)  Same  configuration  as  Figure 

2.22 

2)  The  new  gains  are: 

K_  =  3x10^  N-m/rad/sec 

KAX 

KPxx 

*  3xl07 

N-m/rad 

KD  =  1.05x106  N-m/rad/sec 

ka  y 

KPxy 

=  1 xl U7 

N-m/rad 

KRaz  -  9xl05  N-m/rad/sec 

KPaz 

=  lxl 07 

N-m/rad 

and 

KR37  =  h:106  N/m/sec 

KP37 

*  1 .75x1 08 

N/m 

KR46  =  3x10^  N/m/sec 

KR46 

=  3. 5xl07 

N/m 

/•  .  j _ 

Figure  2.26  Steady  State  Disturbance  Response  of  Tuned-Up  System 
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3.  ON-ORBIT  TEST  AND  SYSTEM  IDENTIFICATION 


3.1  INTRODUCTION 

The  prime  motivation  for  conducting  on-orbit  tests  and  system  identi¬ 
fication  is  the  uncertainties  In  the  finite  dimensional  model  that  approx¬ 
imates  the  dynamics  of  the  large  space  structure.  The  finite  element 
method  is  often  used  to  obtain  this  finite  dimensional  model  which,  from 
past  experience,  is  expected  to  contain  errors  greater  than  20%  in  the 
predicted  modal  frequencies  and  mode  shapes  of  some  of  the  modes.  As  a  conse¬ 
quence,  the  active  controllers  which  are  computed  on  the  basis  of  this 
inaccurate  model  will  not  achieve  their  designed  performance.  Hence  it 
is  desirable  to  obtain  a  more  accurate  model  by  testing  the  actual  space¬ 
craft  on-orbit.  Ground-based  tests  for  large  and  flexible  structures  are 
not  always  feasible;  besides,  the  test  environment  on  the  ground  is  not 
representative  of  on-orbit  conditions.  For  the  above  reasons,  on-orbit 
tests  and  system  Identification  are  important  Issues  in  the  active  con¬ 
trol  technology  for  the  large  space  structures. 

System  identification  is  the  problem  of  determining  the  dynamic 
model  of  an  uncertain  system  based  on  measurements  of  inputs  and  outputs. 

It  Is  a  very  broad  problem  which  even  includes  the  case  where  the  order 
of  the  system  is  unknown.  However  in  this  report,  a  finite  order, 
parametric  description  of  the  dvnamic  models  is  used  where  certain  coeffi¬ 
cients  In  the  dynamic  equations  are  unknown.  A  parame+cr  estimation 
technique  Is  then  used  to  solve  the  system  identification  problem. 

The  concept  of  on-orbit  tests  and  controller  tuning  is  shown  in 
Figure  3-1.  In  this  approach,  LSS  on-orbit  input-output  data  are  trans¬ 
mitted  to  ground  where  the  parameter  estimation  algorithm  is  applied  to 
update  the  LSS  parameters.  The  improved  LSS  model  is  then  used  to  re¬ 
compute  new  control  parameters  in  order  to  tune  the  controller  to  achieve 
the  desired  performance.  A  ground-based  estimator  will  relieve  much  of 
the  computational  burden  from  the  spacecraft  and  allow  more  sophisticated 
parameter  estimation  techniques  to  be  applied. 


CONTROL  PARAMETER  UPDATES 


Figure  3-1.  On-Orbit  System  Identification  and  Controller  Tuning 


3.2  PARAMETER  ESTIMATION  METHODS 


Several  existing  parameter  estimation  methods  are  evaluated  quali¬ 
tatively  for  their  suitabilities  of  performing  on-orbit  parameter 
estimation. 

(A)  Cross-Correlation  Method 

The  error  between  the  measured  output  and  the  computed  output 
is  cross-correlated  with  the  signal  in  the  model  that  Is  an 
input  to  the  parameter  to  be  Identified.  The  output  of  the 
cross-correlation  integral  then  adjusts  the  corresponding 
parameter  value  in  the  system.  The  main  advantage  is  that 
it  can  be  implemented  for  on-line  parameter  estimation.  But 
there  remains  some  difficulties  for  multtple-input-multiple- 
output  systems. 

(B)  Linear  Prediction  Method 

The  parameters  in  an  auto  regression  model  are  selected  so  that 
the  mean  square  prediction  error  is  minimized.  There  are 
efficient  computational  algorithms  to  compute  the  estimates. 
However,  this  method  is  sensitive  to  non-white  measurement 
noise  and  non-zero  initial  states. 

(C)  Transfer  Function  Method 

The  magnitude  square  of  the  LSS  transfer  function  is  obtained 
from  the  ratio  of  the  power  spectral  density  of  measured  out¬ 
put  to  that  of  the  input  signal  through  the  use  of  Fast  Fourier 
Transforms  (FFTs) .  Although  this  method  is  sensitive  to 
non-white  measurement  noise  and  non-zero  initial  states,  it 
is  simple,  easy  to  implement,  and  can  be  used  to  obtain  the 
"first-cut"  solution. 

(D)  Kalman  Filtering  Method 

Parameters  are  treated  as  states  and  the  parameter  estimation 
problem  is  formulated  as  a  non-linear  state  estimation  problem. 
The  extended  Kalman  filtering  algorithm  is  then  applied.  Con¬ 
sequently,  this  algorithm  can  be  used  for  on-line  parameter 
estimation,  but  the  method  is  complex  computationally  and  may 
lead  to  biased  estimates, 


(E)  Maximum  Likelihood  Method 


The  paramettrs  are  chosen  so  that  the  likelihood  functional  Is 
maximized.  When  the  underlying  statistics  is  Gaussian,  the 
maximum  likelihood  method  and  the  least  squares  method  are 
equivalent.  This  method  is  known  to  yield  asymptotically 
unbiased  and  consistent  estimates. 

Among  the  methods  presented  above,  the  maximum  likelihood  method  and 
the  transfer  function  method  are  believed  to  be  most  suitable  for 
on-orbit  parameter  estimation. 

3.3  MAXIMUM  LIKELIHOOD  PARAMETER  ESTIMATION  FOR  LARGE  SPACE  STRUCTURES 

T-/is  section  presents  the  maximum  likelihood  method  [1]  and  its 
applications  in  the  estimation  of  LSS  physical  and  modal  parameters. 

The  effects  cf  data  length  are  also  discussed. 

3.3.1  Maximum  Likelihood  Method 

The  three  major  problems  in  applying  the  maximum  likelihood  method 
to  LSS  are: 

1)  obtaining  an  approximated  finite  dimensional  model  for  the  LSS, 
and  determining  which  parameters  are  to  be  adjusted; 

2)  deriving  an  expression  for  the  likelihood  cost  functional*;  and 

3)  selecting  an  algorithm  to  minimize  the  likelihood  cost  functional 

These  three  problems  and  the  properties  of  the  maximum  likelihood 
estimates  are  discussed  below. 

Modeling  LSS  Dynamics 

The  large  space  structure  is  a  distributed  parameter  system  with 
point  actuators  and  point  sensors.  The  dynamics  of  a  LSS  can  be  described 
completely  by  partial  differential  equations  (PDEs)  in  the  general  form: 


If  an  approximated  model  is  obtained  for  the  LSS,  the  derived  likeli 
hood  cost  functional  will  also  be  an  approximation 
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x(t)  =  AX(t)  +  Bu(t) 
s(t)  =CX(t) 


(3.1a) 

(3.1b) 


where  A  is  the  infinitesimal  generator  of  a  strong  continuous  semigroup 

% 

over  a  separable  Hilbert  space;  X(t)  is  in  D(A);  B  and  C  are  influence 
operators  for  point  actuators  and  point  sensors;  u(t)  is  the  input  vector 
and  s(t)  is  the  ideal  (noise-free)  output  vector.  The  actual  sensor 
measurement  is 

z(t)  =  s ( t )  +  n(t)  (3.2) 

where  n(t)  is  zero  mean  white  Gaussian  noise  with  positive  covariance 
matrix  G. 

Except  for  simple  structures  such  as  thin  plates  and  thin  beams, 
the  exact  solution  of  equation  (3.1)  for  a  complex  LSS  will  be  very 
complicated,  if  not  impossible  to  obtain.  Hence,  one  would  like  to 
approximate  the  dynamics  of  a  LSS  by  a  finite  dimensional  model.  The 
finite  element  method  [2]  is  often  used  to  obtain  such  a  model. 

In  the  finite  element  method,  the  structure  is  discretized  into 
small  basic  elements  connected  at  node  points.  The  motions  of  the  node 
points  can  be  described  by  a  set  of  simultaneous,  second  order  differ¬ 
ential  equations  (with  damping  terms  omitted  for  the  moment), 

Mq(t)  +  Kq(t)  =  Bpu(t)  (3.3a) 

y(t)  =  Cpq(t)  (3.3b) 

where  q(t)  is  the  physical  coordinates  of  the  node  points,  M  is  the  mass 
matrix,  K  is  the  stiffness  matrix,  Bp  and  Cp  are  input  ari'_  output  in¬ 
fluence  matrices,  u(t)  is  input  vector,  and  y(t)  is  the  output  vector 
predicted  from  the  model. 

Let  A  and  <t>  be  the  eigenvalues  matrix  and  the  corresponding  eigen¬ 
vector  matrix  of  the  eigenvalue  problem 

M$a  =  K$ 

Equation  (3.3)  can  be  transformed  into  a  set  of  decoupled  equations 
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where 


(3.4a) 

(3.4b) 


Viscous  damping  can  now  be  modeled  in  terms  of  damping  ratios, 
i  =1,  2,  ...,  n.  Including  damping  terms  and  reordering  the  state 
variables,  one  can  rewrite  equation  (3.4)  as 

x(t)  =  Ax(t)  +  Bu(t)  (3.5a) 

y(t)  =  Cx(t)  (3.5b) 

where 


Equation  (3.5)  is  a  finite  dimensional  approximation  of  the  actual  LSS 
dynamics  described  in  equation  (3.1).  As  mentioned  before,  this  finite 
dimensional  model  is  expected  to  have  errors  in  modal  frequencies, 
damping  ratios,  and  mode  shapes.  The  errors  can  be  minimized  by  adjust¬ 
ing  the  parameters  in  equation  (3.5).  Figure  3-2  shows  two  ways  to  tune 
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Figure  ''-2.  Modal  Parameter  Estimation  and  Physical  Parameter  Estimation 


3-7 


the  model.  The  first  approach  is  to  adjust  the  modal  frequencies, 
damping  ratios,  and  mode  shapes,  called  modal  parameters,  in  equation 
(3.5)  directly.  The  other  approach  is  to  adjust-  the  p?>ameters  in 
matrices  M  and  K,  called  physical  parameters,  in  equation  (3.3). 

Examples  of  physical  parameters  are  Young's  modulus  of  structural  material, 
mass  distributions  on  the  structure,  etc.  Changes  in  M  and  K  will  in¬ 
fluence  the  valies  of  modal  frequencies  and  mode  shapes  via  the  eigen¬ 
value  problem.  Damping  ratios  can  be  adjusted  with  physical  parameters 
at  the  same  time. 

Likel-ihood  Cost  Functional 

Let  the  actual  LSS  dynamics  ani  sensor  measurement  be  described 
by  equations  (3.1)  and  (3.2).  Then  the  likelihood  cost  functional  for 
a  record  of  LSS  input-output  measurements  {u(t),  z(t)}  during  the  time 
interval  [0,T]  is  [3] 

f\j  T 

f  [z(t)  -  s(t)]V[z(t)  -  s(t)]dt  (3.6) 

o 

where  s(t)  is  the  ideal  senso '  measurement  described  in  equation  (3.1). 

The  above  likelihood  cost  functional  which  involves  the  solution  of 
PDE  is  difficult  to  evaluate.  '  w^ver,  it  can  be  approximated  by 

J(e)  =  j  f  [z(t)  -  y(9;t)]TG_1[z{t)-y(e;t)]dt  (3.7) 

0J 

where  y(e,t)  is  the  output  estimated  by  the  finite  dimensional  model 
(equation  (3.5))  which  is  parameterized  as 

xf';t)  =  A(e)  x  (e  ;  t )  B(e)  u  ( t )  (3.8a) 

y'  ,t)  =  c(e)  x (e ;t )  (3.8b) 

and  e  =  [a^,  c^,  •••»  ap^  a  vector  of  unknown  parameters  (modal  or 
physical)  in  matrices  A,  B.  and  C. 

The  maximum  likelihood  estirates  0  is  the  e  which  minimizes  J(e). 

The  -  rameter  estimation  problem  r.ow  becomes  a  non-linear  minimization 
problem .  The  modified  Newton-Raphson  method  can  be  used  to  perform  the 
minimization  numerically. 
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Modified  Newton-Raphson  Method 


Let  denotes  the  a-priori  value  of  the  unknown  pa/amete^  o,  e.g., 
the  nominal  parameter  values.  The  Newton-Raphson  method  is  an  iteractive 
pro».idure  In  which  the  new  estimate  e  .j  is  commuted  from  the  present 
estimate  en: 

Vl  en  -  H_'  <V  <3' 

where  vj(e)  is  the  gradient  of  J(e )  with  respect  toe,  and  H(e)  is  the 
Hessian  matrix.  If  all  the  second  derivatives  in  H(e)  are  omitted,  then 
this  procedure  is  called  the  modified  Newton-Raphson  method. 


9  =  [ai ,  a2* 


-iT 

'»  apJ  . 


then  for  modified  Newton  -  Raphson  method: 


vj(e)  = 


5J(e) 

— 


3 J(e )  =.f  3y(e  ;t)~l  G"1  ]z(t)  -  y(e;t)]  -t 

3a.  I  L  5a.  J  l 


(3. to: 


- T 


where 


ay(o;t) 

9a  - 
1 

is  the  output  sensitivity  with  respect  to  parameter  <*.. 

In  order  for  the  modified  Newton-Raphson  iteration  to  converge,  the 
matrix  H(e)  should  be  positive  definite  in  a  region  containing  the  'rue 
parameter.  This  is  the  identifiabil ity  condition  for  the  maximum  likeli¬ 
hood  method. 

The  output  sensitivity  is  important  in  calculating  the  gradient  and 
the  Hess. an  matrix.  The  computation  of  the  output  sensitivity  for  physical 
and  modal  parameters  will  be  discussed  later. 


Properties  of  Maximum  Likelihood  Estimates 


It  i:  known  that  the  maximum  likelihood  estimates  are  asymptomatically 
unbiased  and  consistent. 


Let  be  the  true  parameter  and  e  be  the  maximum  likelihood  estimate, 
then  the  covariance  matrix  of  the  estimation  error 


where  H  is  the  Fisher  information  matrix.  The  elements  of  H  is  described 
in  equation  13.11).  The  matrix  is  the  Cramer-Rao  lower  bound. 
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(3.13b) 


y( o ; t )  =  Cp  q( 0 ;t) 

0  =  [a1 ,  a2 ,  -  ,  ap]T 


The  output  sensitivity  can  be  obtained  by  differentiating  equation  (3.13) 
with  respect  to  the  unknown  parameter  a.. 


„(e,  .  K(9)  isliil) .  .  | (e;t) 

oClj  dCt^ 


3a  1 


M’1(e)  CK(e)  ^(0;t)  -  V(t)  ] 


(3.14a) 


3y(e;t)  _  r  3q(e;t) 

3a.  "  3tti 


(3.14b) 


For  complex  structures,  the  dimensions  of  equations  (3.13)  and  (3.14) 
could  become  very  large,  and  it  would  bo  numerically  advantageous  to  trans¬ 
form  both  equations  into  sets  of  decoupled  equations  in  the  modal  coordinates. 
Table  3-1  describes  the  transformations.  Equation  (3.13)  is  transformed  into 
modal  coordinates  by  the  modal  matrix.  Except  for  the  input  uart,  equaf'jn 
(3.14)  is  similar  to  equation  (3.13)  hence;  it  can  also  be  transformed  into 
a  set  of  decoupled  equations. 


The  block  diagram  in  Figure  3-3  illustrates  the  use  of  finite  element 
model  in  the  maximum  likelihood  estimation  of  LSS  physical  parameters. 
Initial  or  updated  unknown  physical  parameters  and  other  fixed  parameters 
are  fed  in  a  finite  element  mode  ■  alysis  program  which  generates  new 
system  matrices  for  LSS  dynamics  and  output  sensitivity  computations.  The 
finite  element  program  can  be  a  special  subroutine  written  for  a  particular 
LSS  structure,  or  it  can  be  a  "canned"  program  such  as  NASTRAN  which  will 
be  capable  of  handling  <. nrious  large  and  complicated  structures.  Currently, 
a  special  finite  element  subroutine  is  used.  Also  note  that  when  enough 
mechanical  isolations  are  provided  among  different  parts  of  the  whole 
structure,  the  finitrt  element  p;  jgram  needs  only  to  model  those  parts  whose 
physical  properties  are  of  interest.  This  will  be  demonstrated  later-  in 


»>■***&£■  'il  i  ^  c  ‘-'J-.'.  j- 
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Table  3-1.  Output  Sensitivity  Computation  For  Physical  Parameters 


MhM'K  Nt  AMIKI  MINI  •’ ( t ) 


Maximum  likelihood  with  Newton-Raphson  Iteration  For  Physical 
Parameter  Estimation 


Section  3.6. 


The  damping  ratios,  ?.'s,  are,  strictly  speaking,  modal  parameters, 
but  they  can  be  estimated  with  other  physical  parameters  at  the  same  time. 
This  is  because  the  damping  ratios  are  not  involved  in  the  underlying 
eigenvalue  problem,  and  the  corresponding  output  sensitivities  can  be 
obtained  by  differentiating  the  LSS  modal  equations  with  respect  to 
damping  ratios.  The  damping  ratio  for  each  mode  can  be  treated  as  indivi¬ 
dual  independent  unknown  parameter.  This  will  drastically  increase  the 
number  of  unknown  parameters.  Hence  for  any  practical  applications,  it  is 
desirable  to  divide  damping  ratios  into  small  number  of  groups  so  that 
damping  ratios  in  each  group  are  constraint  to  have  the  same  value.  This 
will  reduce  the  number  of  unknown  parameters.  In  the  subsequent  simulation 
tests,  all  the  damping  ratios  in  the  precision  structure  are  considered  as 
one  unknown  parameter.  The  output  sensitivity  for  this  case  is  computed 
as  follows: 


Then  the  output  sensitivity  for  damping  ratio  can  be  obtained  by  differen¬ 
tiate  equation  (3.5)  with  respect  to  c, 


d  /  3x  (t)\  _  „  3x(t)  +  3A  x(t)  (3.15a) 

dt  \  3c  /  n  3?  a; 


3y(t)  . 
3; 


c 


3xitl 

a; 
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3.3.3  Direct  Estimation  Modal  Parameters 


Accurate  modal  parameters  {modal  frequencies,  mode  shapes)  for  all 
modes  of  interest  can  be  computed  from  the  estimated  physical  parameters 
via  finite-element  modal  analysis  program.  But  the  estimation  of  physical 
parameters  may  be  too  complicated  if  one  only  wants  to  estimate  the  modal 
data  of  a  small  number  of  modes,  or  physical  parameter  estimation  loses 
meaning  because  of  the  complexity  of  the  structure.  The  maximum  likeli¬ 
hood  method  can  be  applied  to  estimate  those  modal  parameters  directly. 

The  output  sensitivities  for  damping  ratios  are  shown  before.  For 
modal  frequencies  and  mode  shapes,  the  output  sensitivity  can  be  computed 
in  a  similar  manner;  i.e., 


d  3x(t)  _  *  3x(t)  .  3A  .  3B 

at  - A 


u(t) 


(3.16a) 


llill  =  C  x  { t ) 

3a  ^  3a  |  3a  1- 


(3.16b) 


where 


a.  is  a  modal  parameter. 

3.3.4  Effect  of  Data  Length 

It  is  known  that  the  accuracies  of  the  estimates  can  be  improved  by 
using  longer  data  measurements.  The  large  space  structure  is  a  very 
lightly  damped  system.  The  corresponding  cost  functional  that  must  be 
minimized  to  obtain  a  good  parameter  estimate  has  a  component  of  the 
sinc-function  type  as  shown  in  Figure  3-4. 

For  long  data  measurements  the  minimum  of  this  function  is  well 
defined,  as  shown  in  curve  (a),  and  the  modified  Newton-Raphson  method 
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will  rapidly  converge,  provided  the  a  priori  (initial)  estimate  eQ  was 
within  the  convex  part  of  the  global  minimum.  Otherwise,  the  algorithm 
will  diverge  to  a  local  minimum  of  th“  quadratic  cost  function.  For  data 
measurements  of  short  length,  the  sine  function  is  broad  as  shown  in 
curve  (b)  and  the  absolute  minimum  is  not  well  defined  for  a  gradient 
search  technique.  Thus,  initially  the  a  priori  parameter  estimates  should 
be  merely  improved  by  processing  relatively  short  pieces  of  measurement 
data.  Thereafter,  the  data  length  should  be  Increased  to  provide  better 
noise  rejection  and  a  "sharp"  convergence  to  more  accurate  parameter 
estimates. 

One  can  expect  that  longer  data  length  will  increase  the  accuracies 
of  the  estimated  parameter,  however,  the  initial  estimates  will  be 
required  to  be  closer  to  the  true  parameter  valves  in  order  for  the  maximum 
likelihood  algorithm  to  converge. 

Simple  Example 

The  effect  of  data  length  can  be  further  illustrated  by  a  simple 
example.  Consider  the  mass-spring  system,  shown  in  Figure  3-5,  with  true 
dynamics 

mQ  x ( t )  +  kQ  x ( t )  =  F(t)  mQ  =  1,  k0  =  1 

and  model 

m  x(t)  +  k  x(t)  =  F(t) 

For  an  impulse  input,  F  =  <S (t) »  the  displacement  of  the  mass  is 

=  A  sin  to  t 
o  o 


x(t)  *  — —  sin  %-  t  =  A  sin  tot 

/IF  * m 


x(t)  = 


1 


/  m  k 
o  o 


sin 


1  m 


Similarly,  the  stimated  response  is  found  to  be 
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T  •  10  SEC  T  •  20  SEC 


3-5.  Cost  Functional  of  Mass-Spring  System  For  D3ta 
Length  T=5,  10  And  20  Seconds  (0.5  <  m  ^  1.5, 
0.5  <  k  1 .5) 


This  cost  functional  can  be  expressed  as 


q(T) 


J  (x(t)  -  x(t))2  dt 


\  »o  *  a2]  -  I  [' 


.2  sin  2“oT  .  ,2  sin  2tuT 
A.  — — =? -  +  A 


0  2“>0T 


2wT 


,  sin  (u  +  oj)T  , 

■  a  a  O  * 


+  2  («0  +  oi)  T 


AA 


sin  (c 


o  (u»0  -  w)T 


The  cost  functional  for  0.5  m  £  1 .5,  0.5  £  k  £  1 .5,  and  T  =  5,  10,  20 
are  shown  in  Figure  3-5. 

!t  can  be  seen  that  as  T  increases  the  cost  functional  becomes 
sharper  around  the  true  values  -  higher  accuracies,  but  one  has  to  start 
at  a  point  closer  to  the  true  values  -  smaller  region  of  convergence. 

3.4  SELECTION  OF  CRITICAL  PARAMETERS 

It  is  known  that  too  many  unknown  parameters  will  decrease  the 
accuracies  of  the  estimated  values  and  even  cause  identifiabil ity  problem 
in  some  cases.  One  need  to  determine  which  parameters  are  critical  and 
must  be  estimated.  All  other  parameters  are  fixed  at  their  nominal  values. 
The  physical  parameters  are  usually  small  in  number  and  hence  pose  less 
difficulties.  Two  criteria  that  can  be  used  to  select  critica1  modal  para¬ 
meters  are  presented  in  this  section. 

3.4.1  Controller  Performance  Sensitivity  Criterion 

Consider  the  modal  representation  of  a  LSS  dynamic  model : 


where  Xc  is  the  controlled  mode  and  Xr  is  the  residual  mode. 

The  actuator  input  signal  is  u(t),  and  the  sensor  measurement  is 

y(t).  For  the  following  analysis,  it  is  assumed  that  the  actuators  and 

sensors  are  collocated,  hence  bI  =  C  and  B*'"  =  C  . 

i*  c  r  r 

The  feedback  controller  selected  for  the  sensitivity  study  is  the 
one  that  satisfies  both  positivity  design  and  LQG  design.  This  particular 
design  may  not  meet  the  controller  performance  requirement,  but  it  pro¬ 
vides  the  generic  structure  of  both  positivity  and  LQG  designs. 


;c  ■  <*c  -  Gc  Cc)  ;c  ♦  Gc  y 
«  -  •  *c 


where 


G  *  B. 

c  c 

*T  * 

Kc  *  GCH 
P  (Ac  -  Gc  Cc)  ♦  (A 

Q  ■  Olag  [Qj] 

0 

0  “j. 


GcCc)  P 


-  [0  ♦  26c  B'] 


'j 


The  symbol  is  used  to  denote  those  quantities  that  are  function  of 
unknown  parameters. 

The  controller  equations  and  the  LSS  dynamics  can  be  combined  to  form  an 
augmented  system. 


A  X 


X  3  AX  X(o)  3  X„ 

o 


The  controller  performance  index  (not  the  performance  Index  for  controller 
design)  is  taken  to  be 


J 


R  X(t)  dt 


Rf  R  eAt  dt  Xq 


Where  the  matrix  R  can  be  chosen  so  that  RX  represents  the  quantity  to  be 
minimized,  e.g.,  RX  could  be  1 ine-of-sight  error.  The  performance  index, 
J,  is  a  function  of  both  A  and  XQ.  A  usual  mathematical  procedure  to 
eliminate  the  dependence  on  the  initial  states  xQ  is  to  average  J  over  the 
n-dinensional  unit  sphere  of  initial  states.  This  leads  to  the  average 

A 

performance  index  0 


J  *  n  E  {J}  3  Trace 


|j  eftTt  RT  R  eAt  dt  =  .'race  [  L  ] 


where  L  is  the  solution  of  a  Liapunov  equation 


AT  L  +  L  A  +  R  R  =  0 


Let  a.,  «2»  - »  ap  be  moda1  parameters,  then  the  sensitivity  of  per¬ 


formance  index  0  to  parameter  a.  is 


Trace 


[’“J 


i  3  1,2 . p 
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where  31  satisfies 
3a . 


at  ik_  +  i-L.  A  + 

da^ 


U! 

hi 


L  +  L 


Note  that  for  physical  parameters  the  computation  o*  sensitivity  is  more 
complicated,  due  to  the  fact  that  the  controller  is  designed  through  the 
use  of  modal  representation.  For  physical  parameters  8-j ,  Bj.  ....  8^, 
the  sensitivity  of  performance  index  is 


L  -  y*  3J  3a1 
1  J  mmm4  3a  t  3  B  1 
J  1  0 


3  t 

where  is  sensitivity  to  modal  parameters,  and  ~  involves  eigenvalue 
da  •  dpj 

problems  but  can  be  approximated  by  using  Jacobi's  formula. 

The  performance  sensitivity  of  each  parameter  is  itself  an  indica¬ 
tion  of  importance,  however  one  can  not  compare  the  performance  sensitivity 
of  mode  frequency  with  the  performance  sensitivity  of  damping  ratio. 

Let  AJ  be  the  total  allowed  degradation  in  the  performance  index 


^  Aal  +  f^Aa2  +"*  +  l^Aap 


where  Aa.j  is  the  change  in  parameter  a^ 


The  degradation  of  AJ  will  not  be  exceeded  if  Aa.'s  are  such  that 


lA.,!  < 

1  3J  I 


i  =  1 ,2, . . .  ,p 


i  parameter  values  are 

|AJ, 
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CONST. 

\lLa 

|3o1  i 

a. 
9a^  i 

Therefore,  in  order  to  assure  controller  performance  robustness,  those 

IActi 

values  (or  equivalently  high vaiues) 

should  be  selected  for  parameter  estimation.  For  instance,  if  one  wants 

to  estimate  m  parameters,  then  the  first  m  parameters  with  highest 

9 J  g.  values  should  be  estimated. 

3a^ 


Example  (DRAPER  Tetrahedral  Truss  Structure) 

This  parameter  sensitivity  analysis  technique  cf  selecting  critical 
system  parameters  is  applied  tc  the  Diaper  tetrahedral  truss  structure 
(See  Section  4-2) 

Controlled  Modes:  Modes  1-8 

Residual  Modes:  Modes  9-12 


bl,l  bl  ,2 


b2 ,1  b2,2 


-2<i  “i 


-  *  T 

c  *  B1 
c  c 


ft 


Performance  Index:  RX 


The  modal  parameters  are  the  modal  frequencies  (t». ),  the  mode 
shapes  (b.,.:)  and  the  damping  ratios  (  ?.)  in  the  controlled  modes. 

'  w  1 


Table  3-2  lists  the  first  21  parameters  (out  of  64)  according  to  their 


3  J 


a  • 

3a-  1 


values.  Large  values  in  the  last  column  means  significant 


influence  on  performance. 


3.4.2  Residual  Power  Spectral  Density  Criterion 

In  previous  criterion,  the  critical  parameters  are  selected  indepen¬ 
dent  of  the  actual  input-output  measurements  of  the  LSS.  Some  parameters 
may  be  very  close  to  their  true  values  already  and  need  not  be  estimated. 
Another  criterion  of  selecting  critical  parameters  is  to  examine  the  Power 
Spectral  Density  (PSD)  plots  of  the  residuals  (output  errors).  Taos*  modes 
whose  frequencies  are  closest  to  the  spikes  in  the  residual  PSDs  are 
chosen  for  parameter  estimation.  One  mode  or  several  modes  can  be  estimated 
at  the  same  time,  and  update  the  modal  data.  This  process  can  be  repeated 
until  desired  error  levels  are  achieved. 


The  application  of  this  criterion  is  illustrated  later  on  the  DRAPER 
Model  #2. 


3.5  DESIGN  OF  ON-ORBIT  TEST  SIGNALS 

In  general,  the  input  signal  to  a  LSS  should  be  "sufficiently  rich" 
so  that  the  modes  of  interests  are  excited.  A  more  systematic  approach  is 
to  examine  the  resulting  Fisher  information  matrix  H  (equations  (3-11)  and 
(3.12)).  The  inverse  of  H  is  an  accuracy  measure  on  the  estimated  para¬ 
meters.  It  would  be  beneficial  to  find  an  input  signal  which  optimizes 
functionals  of  H  subject  to  some  energy  constraint.  Examples  of  functionals 


s~ 

*  i 


n 


M 
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are:  (Det(H),  Trace  (H),  and  Trace  This  is  the  optimal  input 

signal  design  problem. 

The  optimal  input  signal  can  assure  the  identifiabil ity  condition 
and  improve  the  estimation  accuracies.  In  other  words,  the  same  levels  of 
estimation  accuracies  can  be  achieved  with  shorter  LSS  input-output 
measurement  data.  Hence,  the  time  required  in  on-orbit  testing  and  para¬ 
meter  estimation  can  be  reduced. 

But  optimal  input  design,  in  general,  is  a  difficult  problem  to  solve. 
The  results  depend  on  the  unknown  parameters,  and  may  violate  the  space¬ 
craft  mission  requirement,  structure  load  constraint,  or  linearity  assump¬ 
tion.  An  alternative  approach  to  input  design  problem  can  be  as  follows. 
Several  input  signals  which  satisfy  the  mission  requirement  and  the 
structure  load  constraint  are  selected  as  candidates.  The  one  with  highest 
Det(H)  value  is  used  to  excite  the  LSS. 

3.6  APPLICATIONS  TO  DRAPER  MODEL  #2 

Computer  simulation  tests  were  conducted  to  verify  the  performance 
of  the  maximum  likelihood  parameter  estimator  as  it  applied  to  the  Draper 
Example  Structure  Model  #2  (see  Appendix  A).  The  results  of  estimating 
physical  parameters  and  modal  parameters  are  presented. 

3.6.1  Draper  Model  # 2 


The  Draper  example  structure  Model  *2  and  its  actuator  and  sensor 
locations  are  shown  in  Figure  3-6.  The  actuator  is  of  force  type,  and  it 
can  only  act  in  the  truss  axial  direction.  The  sensor  measures  the  relat¬ 
ive  position  displacement  of  two  nodes  in  the  direction  of  the  line 
connecting  them.  Two  collocated  actuator./sensor  pairs  are  placed  on  truss 
members  connecting  nodes  27  to  30,  and  9  to  10.  Two  additional  sensors  are 
from  nodes  2S  to  32  and  28  to  35  which  has  no  connecting  truss  members.  The 
actuator  and  sensor  from  node  27  to  30  is  denoted  as  "actuator  27-30"  and 
"sensor  27-30",  respectively. 
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Note  that  the  actuator/sensor  locations  differ  from  that  of  the 
active  control  (Section  2).  This  is  because  the  control  studies  and 
identification  studies  are  conducted  independently  at  the  same  time,  in 
practice,  one  should  attempt  to  perform  parameter  estimation  using  the 
same  set  of  actuators  and  sensors  that  are  selected  for  control.  If  the 
actuators  anc  sensors  for  control  are  not  sufficient  for  identification, 
additional  actuator/sensor  should  be  placed. 

The  perturbed  Model  #2  is  taken  to  be  the  true  system.  The  input 
signal  for  actuator  27-30  is  shaped  pulse  (100Hz  bandwidth)  of  magnitude 
1  Newton,  duration  .02  second,  and  period  .25  second.  Same  input  sig.ial 
with  .02  second  time  delay  is  applied  to  actuator  9-10.  The  simulated 
sensor  measurements  are  generated  by  using  the  perturbed  Model  #2  which 

n 

contains  84  modes.  A  small  amount  (o  =  3  x  10  m)  of  independent  white 
Gaussian  noise  is  added  to  the  simulated  sensor  measurements.  The  highest 
modal  frequency  is  at  247  Hz,  but  the  modes  of  interests  have  modal 
frequencies  below  GO  Hz.  The  sampling  rate  is  taken  to  be  1000  Hz  Which  is 
fast  enough  so  that  the  continuous  time  formulation  for  the  maximum  likeli¬ 
hood  estimator  is  essentially  valid  for  this  set  of  data. 


3.6.2  Physical  Parameter  Estimation 
Problem  Statement 


The  objective  is  to  extract  physical  parameters  from  the  simulated 
data.  The  unknown  parameters  are  the  Young's  modulus,  E,  of  all  truss 
members  in  the  precision  section  and  the  mass  distributions  at  upper  and 
lower  support  structures.  The  initial  and  true  parameter  values  are 
listed  in  Table  3-3.  In  order  to  reduce  the  number  of  unknown  parameters 
the  masses  at  some  nodes  are  constraint  to  have  identical  weights.  The 
shearing  modulus  of  elasticity,  G,  is  also  constraint  so  that  G=0.3E. 

It  is  assumed  that  the  estimator  has  correct  information  on  the  values 
of  A,  I,  and  J  for  all  truss  members. 

The  damping  ratio  was  fixed  at  its  true  value  in  the  intial  tests 
and  was  estimated  alone  with  other  unknown  parametc  in  later  test  runs. 


Simplified  Structural  Model  for  Physical  Parameter  Estimation 


Since  precision  section  and  equipment  section  are  connected  by 
isolators  and  only  the  structural  properties  of  the  precision  section  are 
of  interest,  the  finite  element  model  for  the  e-timator  can  be  simplified 
to  the  model  shown  in  figure  3-7.  The  equipment  section  and  solar  panels 
are  omitted.  The  precision  section  is  mounted  on  isolators  to  inertial ly 
fixed  rigid  support.  The  mode  frequencies  of  this  simplified  model  and 
the  complete  model  are  also  shown  in  the  figure.  All  the  mode  frequencies 
of  the  simplified  model  (except  6  isolator  modes)  can  be  found  in  the 
mode  frequencies  of  the  complete  model  (denoted  by  "a").  In  fact,  this 
simplified  model  keeps  all  the  modes  associated  with  the  precision  section 
intact  while  reducing  the  number  of  modes  from  84  to  54. 

Simulation  Results 

A  total  of  14  test  runs  were  made.  Runs  1  through  11  used  .1 
second  of  simulated  data.  Runs  12  through  14  used  .5  second  of  simulated 
data.  Damping  ratio  was  fixed  (not  estimated)  for  runs  1  through  10, 
and  was  estimated  alone  with  other  unknown  parameter  for  the  last  4  runs. 
The  results  are  summarized  in  Table  3-4.  Each  individual  run  is  high¬ 
lighted  below. 

(A)  Runs  1-5 

Run  1  used  .1  second  of  data,  "'he  estimation  algorithm  starting 
from  the  standard  initial  values  converged  in  6  iteration.  Figure  3-8 
compared  the  time  histories  of  the  simulated  measurements  with  the  estima¬ 
ted  measurements  generated  by  the  estimator  using  standard  initial  para¬ 
meter  values.  Similar  comparisons  were  shown  in  Figure  3-9  in  which  the 
estimated  measurements  were  obtained  from  the  final  estimates  of  Run  1. 

It  appears  that  the  estimator  has  made  significant  improvements  in  the 
estimated  sensor  measurements.  The  residuals  (differences  of  simulated 
and  estimated  measurements)  of  Figure  3-°  are  shown  in  Figure  3-10.  It 
can  be  seen  that  the  residuals  consist  primarily  of  sensor  noise. 
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Figure  3-10.  Residuals  Of  Run  1 
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Figure  3-12.  Corresponding  Modal  Frequencies  For  Each  Iteration  in  Run  1 


The  convergence  of  the  estimated  parameter  In  each  iteration  is 
plotted  in  Figure  3-11.  The  numerical  values  are  listed  in  Table  3-5. 

The  estimated  parameter  converge  after  6  iterations. 

Although  the  mode  frequences  and  mode  shapes  are  not  estimated 
directly,  they  can  be  computed  from  the  estimated  physical  parameters  via 
modal  analysis  program.  The  corresponding  mode  frequencies  for  each 
iteration  are  shown  in  Figure  3-12. 

To  test  the  ability  of  the  identification  algorithm  to  converge 
from  other  intial  parameter  values,  RUNS  2,  3,  4,  and  5  are  made  on  the 
same  simulated  data.  RUNS  2,  3,  and  4  converge  in  6  to  8  iterations. 

RUN  5  which  is  the  worst  case  (up  to  5056  error  in  Initial  parameter  values) 
converged  in  12  iterations. 

(S)  Runs  6,  7 

In  Run  6,  12  highest  frequency  modes  are  truncated.  This  has  very 
little  effect  on  the  estimated  values.  However  in  RUN  7,  with  18  highest 
frequency  modes  truncated,  there  are  noticeable  degradations  in  the  esti¬ 
mated  parameter  values.  It  is  possible  to  truncate  more  modes,  but  it 
has  to  be  done  selectively,  l.e.  dominant  modes  should  be  retained. 

(C)  Runs  8-10 

In  RUN  8  all  damping  ratios  are  to  be  .05%  while  true  values 
are  .156.  There  are  very  little  changes  in  the  estimates,  and  this  indicates 
that  the  damping  rations  will  be  difficult  to  estimate  from  this  set  of 
data. 


In  RUN  9,  the  estimator  has  Incorrect  values  for  the  masses  In  the 
metering  truss.  RUN  10  simulates  when  sensor  28-35  fails,  and  the  measure¬ 
ments  from  it  are  not  used  in  the  estimator.  Both  runs  show  some  degrada¬ 
tions  in  the  estimated  parameter  values. 


(D)  Run  11 


In  RUN  II,  the  damping  ratio  ;  (for  all  structural  modes)  Is 
estimated,  in  addition  to  other  physical  parameters,  using  0.1  second  of 
data. 


It  is  assumed  that  the  value  of  ?  is  known  to  lie  between  0.01% 
and  1.00%.  The  estimator  converges  in  9  Iterations.  The  damping  ratio 
is  the  last  parameter  to  converge  and  its  final  estimate  is  .075%  (true 
value  =.100%).  But  according  to  the  C-R  Bound  (standard  deviation  =.05%) 
this  final  estimate  is  a  very  reasonable  value  for  0.1  second  of  data. 

To  obtain  better  estimates,  longer  data  are  to  be  used. 

(E)  Runs  12-14 


To  obtain  more  accurate  estimates,  the  maximum  likelihood  estimator 
is  applied  to  0.5  second  of  simulated  data.  The  first  0.1  second  of  this 
data  is  identical  to  the  data  used  in  previous  test  cases. 

RUN  1?  has  the  same  initial  parameter  values  as  RUN  11,  but  after 
9  iterations,  the  estimator  shows  no  signs  of  convergence.  The  divergence 
is  caused  by  too  large  errors  in  the  initial  parameter  values  relative  to 
data  length  (discussed  in  Section  3.3.4).  Hence,  in  RUN  13,  the  final 
estimate  of  RUN  11  are  used  as  the  initial  estimates.  The  estimator 
converges  in  6  iterations.  All  the  final  estimates  are  very  close  to  the 
true  values,  especially  for  the  damping  ratio.  The  time  histories  of 
simulated  and  estimated  measurements  are  compared  in  Figure  3-13. 

The  initial  parameter  values  for  RUN  13  are  very  close  to  the  true 
values  already.  RUN  14  represents  a  case  between  two  previous  runs  as  far 
as  initial  parameter  errors  are  concerned.  The  estimator  converges  in  9 
iterations.  The  final  estimates  are  the  same  as  RUN  13. 
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Figure  3-14.  Comparisons  Of  Simulated  Measurements  With  Measurements 
’  Predicted  By  3-Mode  Model  Using  Nominal  Modal  Data. 
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Figure  3-15.  Residual  Power  Spectral  Density  Functions  For  Nominal  Hodal  Data. 


3.6.3  Direct  Estimation  of  Modal  Parameters 

The  objective  is  to  improve  the  nominal  modal  data  by  estimating 
the  modal  parameter  of  few  critical  modes  which  are  selected  via  the 

residual  power  spectral  density  (PSD)  function  criterion  described  in 

Section  3.4.2. 

Simulation  Results 

The  estimator  has  a  dynamic  model  containing  30  dominant  modes: 
modes  17,  21,  28,  30,  34-49,  51-54,  and  56-61.  The  simulated  measurements 
(.5  second)  are  the  same  as  before.  Initially,  all  30  modes  in  the  estima¬ 
tor  are  set  to  their  nominal  values.  The  estimated  measurements  from 
nominal  values  are  compared  to  the  simulated  measurements.  The  time 
histories  are  shown  in  Figure  3-14.  The  residuals  PSD  functions  are 
shown  in  Figure  3-15.  From  the  residual  PSD  plots,  all  sensors,  except 
sensor  9-10,  show  large  errors  at  frequencies  of  100  rad/sec  (mode  41)  and 
210  rad/sec  (mode  58). 

The  modal  parameters  of  mode  41,  except  C£  ^ ,  are  estimated.  The 
new  residual  PSD  plots  are  examined,  and  mode  58  is  chosen  for  further 
parameter  estimation.  After  modes  41  and  58  are  identified,  the  time 
histories  comparisons,  and  the  residual  PSD  plots  are  shown  in  Figure  3-16, 
and  3-17,  respectively.  It  appears  that  the  maximum  likelihood  estimator 
reduces  the  residuals  by  a  factor  of  10  in  power  (3  in  magnitude). 

The  numerical  values  of  final  estimates  for  modes  41  and  58  are 
listed  in  Table  3-6.  The  estimates  are  not  as  close  to  their  true  values 
as  seen  in  previous  cases  where  physical  parameters  ire  estimated.  This 
is  because  all  other  parameters  are  at  their  erroneous  nominal  values,  and 
exact  convergence  cannot  be  expected.  Nevertheless,  a  very  significant 
improvement  in  the  knowledge  of  critical  modes  has  been  achieved. 

Damping  ratio  cannot  be  estimated  to  the  desired  accuracy  from 
this  set  of  data.  However,  a  run  is  made  where  all  modes  are  artificially 
set  to  their  true  values,  except  the  ones  to  be  estimated,  all  parameters 
including  damping  ratios  converge  to  the  true  values. 
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SYSTEM  DESIGN  CRITERIA 


The  operations  for  proposed  on-orbit  parameter  identification  are: 
collecting  LSS  data  on-orbit,  transmitting  data  to  ground,  and  estimating 
data  on-ground.  Since  estimations  are  done  on-ground,  there  will  be  only 
moderate  software/hardware  requirements  for  the  spacecraft.  The  areas 
pertaining  to  system  design  are  in  control  system  interfacing,  actuator/ 
sensor  placements,  actuator/sensor  dynamics  and  data  gathering  and 
handling: 

•  Identification  Phase 

It  should  be  determined  whether  input-output  data  obtained 
from  LSS  normal  maneuvers  are  sufficient  for  parameter  identi¬ 
fication.  If  identification  should  be  included  as  a  phase  of 
LSS  operations,  it  should  be  decided  that  the  control  loop 
be  open,  closed,  or  closed  with  reduced  gain. 

•  Actuator/Sensor  Placement 

Check  whether  the  number  of  actuators/sensors  for  control 
system  are  sufficient  for  parameter  identification.  If  not, 
where  and  what  type  of  additional  actuators/sensors  should 
be  placed. 

0  Actuator/Sensor  Bandwidth 

Actuator/sensor  bandwidth  is  preferred  to  be  much  higher  than 
the  i-^irtant  mode  frequencies.  If  bandwidth  and  mode  fre¬ 
quencies  are  in  the  same  order  of  magnitude,  the  actuator/ 
sensor  transfer  functions  should  be  known  to  high  degrees  of 
accuracies.  Sensor  noise  level  should  be  known  within  an 
order  of  magnitude. 
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Input  Signal 


Input,  signals  should  be  "sufficiently  rich"  (e.g.,  high  Det(H) 
value)  and  still  within  safety  constraint. 

•  Data  Handling 

Tradeoffs  between  data  length  and  estimation  accuracies 
should  be  studied.  Other  topics  such  as  the  quantization 
level,  sampling  rate,  data  storage  space,  and  synchronization 
among  the  signals  of  actuators  and  sensors  should  also  be 
determined. 

3.8  CONCLUSIONS 

The  maximum  likelihood  method  presented  in  this  section  has  shown 
significant  premise  for  estimating  LSS  parameters  from  input-output  measure¬ 
ment  data  obtained  on-orbit.  Simulation  studies  were  conducted  for  both 
physical  parameter  and  modal  parameter  estimation.  Excellent  results  for  estimating 
physical  parameters  were  obtained.  Significant  improvements  on  nominal 
modal  data  were  also  achieved  by  estimating  the  modal  parameters  directly. 

Modal  parameters  require  a  relatively  simple  algorithm  for  estimation,  bu‘: 
they  are  large  in  number  and  become  "over-parameterized"  easily.  On  the 
other  hand,  the  physical  parameters  are  usually  small  in  number,  but  the 
identification  algorithm  (which  involves  finite  element  modeling)  is  more 
compl icated. 

Another  important  subject  studied  was  the  effects  of  measurement 
data  length.  The  accuracies  of  estimated  parameters  can  be  improved  by 
using  longer  data  measurements.  But  it  is  known  in  theory,  and  is  also 
observed  in  the  simulation  tests  (RUNS  11,  12,  and  13),  that  for  longer 
data,  the  initial  parameter  estimates  have  to  be  closer  to  their  true 
values  in  order  for  the  maximum  liNelihood  estimator  to  converge. 
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COMPARISON  OF  DIRECT  ADAPTIVE  CONTROL  WITH  FIXED  GAIN  CONTROL 


The  study  of  direct  adaptive  control  is  motivated  by  its  ability  to 
account  for  system  changes  and  modelling  errors.  These  errors  arise  in  LSS 
control  through  modelling  inaccuracies,  component  failures,  depletion  of 
consumables,  changes  in  the  working  environment  (e.g.,  thermal),  and  through 
reorientation  of  possibly  large  articulated  payloads.  In  most  cases,  the 
on-orbit  identification  procedure  described  in  Section  3  can  be  combined  with 
standard,  fixed  gain,  control  techniques*  to  solve  these  problems.  Direct 
adaptive  controls  may  provide,  however,  an  on-board  autonomous  alternative 
which,  under  some  circumstances  (e.g.,  extended  periods  of  no  communication 
with  ground  links),  is  more  desirable. 

In  this  section  direct  adaptive  control  will  be  compared  with  two 
fixed-gain  control  approaches.  Positivity  constrained,  direct  adaptive 
control  theory  will  be  summarized  first.  Then  fixed-gain  control,  direct  adaptive 
controls  and  identification  followed  by  fixed-gain  control  will  be  compared. 

The  approach  will  be  to  use  the  Draper  Tetrahedron  as  representative  of  a 
large  space  structure;  design  a  control  system  for  each  type  of  control  such 
that  equivalent  performance  under  nominal  conditions  is  obtained;  make  perfor¬ 
mance  comparisons  using  a  perturbed  model  . 

4.1  POSITIVITY  CONSTRAINED  ADAPTIVE  CONTROL  THEORY 

This  sectijn  will  outline  the  theory  behind  positivity  constrained 
adaptive  control.  First,  the  required  extension  of  the  positivity  theorem 
used  in  Chapter  2  will  be  described.  Second,  since  digital  control  will  be 
the  rule  rather  than  the  exception  in  future  applications,  the  discrete 
time  analogs  to  the  concepts  used  in  Chapter  2  will  be  summarized.  Third, 
it  will  be  shown  how  the  adaptive  control  used  here  evolves  naturally  from 
well  established  principles.  Finally,  it  will  be  shown  how  the  adaptive 
controller  is  mechanized  so  that  it  can  be  implemented  with  current  micro¬ 
processor  technology. 


*  This  design  approach  is  generally  termed  "indirect"  adaptive  control. 
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4.1.1  Positivity  Theory  For  Non-Linear  Systems 


The  adaptive  controls  to  be  considered  here  are  highly  nonlinear 
and  time-varying  and  therefore  require  a  more  generalized  approach  to  esta¬ 
blish  stability.  As  alluded  to  in  Chapter  2,  positivity*  is  a  general 
concept  not  limited  to  linear  time-invariant  systems,  and  is  therefore  very 
able  of  treating  the  problem  here. 

The  following  definitions  will  be  required  to  state  the  general 
positivity  theorem: 

Definition  1:  Truncated  Function: 

f  (t)  =  0<  t  <  T 

(0  ,  t  >  T 

Definition  2:  Extended  Hilbert  Space^g:  The  space  of  all  functions  f(t) 
such  that  fy(t)  belongs  to  the  Hilbert  space  i.e., 

f  ■=?  fT  tih 

Definition  3:  Positive  operator:  An  operator  with  a  domain  and  range  in 
an  extended  Hilbert  space  is  positive/strictly  positive 
if  for  fa 

<  fj,  Hfl>l  0  (>0) 

Positive  operators  are  useful  because  they  have  the  following  pro¬ 
perties: 

1)  The  inverse  of  a  positive  operator  is  positive. 

2)  The  sum  of  positive  operators  is  positive. 

3)  Positivity  Theorem  [1,  12]:  The  negative  feedback  systems 
shown  in  Figures  3-la  and  3-lb  are  bounded-input-bounded 

*  Recently  Willems  [9,10]  and  others  [11]  have  made  further  extensions  to 
the  theory  by  using  the  more  general  concept  of  energy  dissipation. 
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Figure  4 


VECTOR  VALUED  SYSTEM  "S" 


lb  Unity  Negative* Feedback  System 
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output  stable  for  all  Inputs  in^p  if  both  G  and  H  are  positive  operators 
and  at  least  one  of  them  is  also  strictly  positive. 

Like  the  positivity  theorem  in  Section  2,  the  advantage  of  this  test 
is  that  it  permits  one  to  Infer  closed-loop  stability  from  the  Individual 
properties  of  the  controller  and  plant.  Also,  like  in  Chapter  2,  when  the 
theorem  cannot  be  applied  directly,  the  same  embedding  techniques  discussed 
earlier  can  be  used  to  establish  stability. 

Finally,  we  note  in  passing,  that  the  term  "Positive  Real"  refers 
to  a  positive  operator  that  is  also  linear  time-invariant,  and,  the  term 
"Hyperstability"  refers  to  the  positivity  theorem  when  at  least  one  of  the 
subsystems  G  or  H  is  positive  real,  while  the  other  may  be  nonlinear  and/or 
time-varying. 

4.1.2  Modifications  Imposed  By  Digital  Positivity  Designs 

Any  control  system  to  be  considered  for  LSS  control  will  likely  be 
digital.  This  is  particularly  true  for  adaptive  control  designs  where, 
because  of  the  algorithmic  complexity,  the  cost  and  flexibility  advantages 
offered  by  microprocessors  are  prevalent.  The  controller  comparisons  made 
here,  therefore  assume  digital  adaptive  and  fixed-gain  controls. 

The  generality  offered  by  positivity  theory  permits  stability 
assurance  of  digital  designs  without  additional  technical  difficulties.* 

For  example,  Anderson's  discrete  time  analog  to  the  continuous  time  domain 
positivity  test  Is: 

Let  A,  B,  C,  D  be  a  minimal  realization  of  G  (z). 

i.e.,  define 

G(z)  =»  D  +  C  (z  I -A)’1  B 


*  The  inner  product  ,  H  f^  merely  changes  from 

f  °f*(t)  *  [Hf { t ) ]  dt  to  jj  f*(iT)  [Hf(1T)]  where  nT«To. 
Jo  1«0 
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where 


then  G(z)  is  positive  real  if  and  only  if  there  exists  a  symmetric 
L  and  W  such  that 

i)  AT  PA  -  P  =  -LLT 

ii)  ATPB  =  CT  -LW 

iii)  W^W  =  o  +  DT  -  BTPB 


Similarly,  the  analogous  frequency  domain  positivity  index  for 
sampled  data  systems  is: 

*{■}  =  >  min  |  \  CG(z)  +  G*(z)]|  ,  z*  ejuT,  «e[o,  Zir) 

where 

£ (•*))  is  always  >0  implies  G(z)  is  strictly  positive  real 

o(u)  is  always  >0  implies  Gfz)  is  positive  real 

5(x)  can  be  <0  implies  G(z)  is  not  positive. 

The  stability  of  sampled  data  systems  is  therefore  handled  by  using  the 
stated  definitions  in  the  positivity  theorem.  Design,  including  embedding, 
is  handled  with  methods  directly  analogous  to  the  methods  of  Chapter  2. 

Unlike  in  the  continuous  case,  discrete  time  systems  almost  always 
require  embedding.  This  occurs  because  sampling  plus  the  zero-order-hold 
introduces  sufficient  phase  lag  to  destroy  positivity,  even  when  ideal 
actuators  and  rate  sensors  are  assumed.  It  is  noteworthy,  however,  that 
once  embedded,  aliasing  does  not  affect  stabil ity  as  long  as  the  system  is 
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G  REPRESENTS  LSS 


observable  and  controllable. 


4.1.3  Evolution  of  Positivity  Design  To  Positivity  Constrained 

Adaptive  Control 

The  outgrowth  of  the  adaptive  design  from  the  positive  design  is 
straightforward.  Figures  2.2a-d  show  the  evolution  by  using  the  positivity 
theorem  and  established  linear  n.odel  reference  control  principles. 

Figure  4.2a  shot's  the  standard  feedback  connection  of  two  a  .s 
G  and  .  G  is  assumed  to  be  the  LSS  and  is  assumed  to  be  the  oiler. 
From  positivity  theory,  we  know  that  if  G  and  are  (strictly)  positive, 
then  the  closed  loop  system  is  stable*. 

In  Figure  4.2b,  an  pre-filter  is  appended.  It  is  well  known 
that  this  so  called  "two-degree-of- freedom  structure"  is  the  most  general 
feeback  controls  configuration  available  and  that  all  other  control  confi¬ 
gurations  can  be  equivalenced  to  it  with  the  proper  choice  of  and  H 2 
[13].  Note,  the  stability  of  this  system  is  guaranteed  if  is  stable 
and  G  and  H1  satisfy  the  positivity  theorem  conditions. 

Figure  4.2c  adds  a  parallel  reference  model  to  the  two-degree-of 
freedom  structure.  It  is  known  that  if  G  and  the  reference  model  are 
linear  and  pre- specified,**  operators  and  exist  which  will  make 
the  error  e  identically  zero  [14].  Note  that  the  reference  model  does  not 
affect  stability  and  in  fact  has  no  function  other  than  to  specify  the 
desired  behavior. 

In  practice,  G  is  highly  uncertain  and  of  unknown  order,  so  and 
cannot  be  pre-specified  such  that  e  is  zero  (or  even  at  a  minimum) . 

Suppose  H-j  and  ^  are  continuously  chosen  to  minimize  e  but  subject  to 


★ 


It  is  understood  that  if  G  is  not  positive,  then  embedding  is  used. 
G  must  be  of  known  order  and  minimum  phase 


positivity  conditions  based  on  a  priori  knowledge  (albeit,  crudely)  of  G. 

In  this  case,  performance  could  be  improved  on-line  without  sacrificing 
the  8H-.  -  loop  stability.*  The  idea  expressed  is  graphically  shown  in 
Figure  4. 2d  and  will  be  called  positivity  constrained  adaptive  control. 

Note  that  because  of  the  generality  of  the  concepts  used,  the  system  can  be 
continuous  or  digital. 


Adaptive  Controller  Mechanization 


The  general  structure  for  the  positivity  constrained  controller  has 
been  shown  in  Figure  4-2d.  The  actual  mechanization  of  the  positivity 
constraints  and  the  adaptive  control  algorithm  will  now  be  described  and 
the  implied  computational  requirements  summarized. 


Testing  The  Positivity  Of  The  Compensators  and  One  of  the 

reasons  for  using  adaptive  control  is  that  the  LSS  dynamics,  G,is  uncertain 
and  possibly  time-varying.  Therefore,  the  compensators  and  (referring 
to  Figure  4. 2d)  are  subject  to  change  on-orbit  and  the  adaptive  algorithm 
must  include  checking  the  positivity  of  H-j  and  H 2  to  assure  stability.  This, 
in  general  can  be  done  using  one  of  Anderson's  positivity  tests.  However, 
it  would  be  a  very  complex  task  and  could  hardly  be  performed  on-line  in  real 
time,  unless  the  configurations  of  and  are  somehow  restricted.  The 
idea  here  is  to  let  the  compensators  he  diagonal.  The  H11  element  can  then 
be  put  in  the  form** 


That  .s,  each  diagonal  element  of  the  compensators  consist  of  a  parallel 
bank  of  filters  and  a  frequency  independent  gain  K. ..  Since  the  sum  of 
positive  operators  is  positive,  then  if  each  of  the  terms  (parallel  filters) 
in  H11  is  strictly  positive,  then  H11  is  strictly  positive  and  so  is  H. 

*  Guarantee  of  a  bounded  parameter  update  sequence  is  required  to 
establish  adaptive  loop  stability. 

**  D  is  the  derivative  operator.  For  digital  systems,  replace  D  by  the 
delay  operator  2. 
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Checking  the  positivity  of  H  is  now  straightforward,  requiring, 
in  the  continuous  case,  that  each  second  order  section 


a-jD  +  a2 


D  +  b}D  +b2 


be  constrained  by 


i)  a-j,  a2,  b^ ,  b2^0 
fi)  a,b,>a2 


Analogously,  in  the  discrete  time  case,  each  second  order  section 


a-j  2  +  a2Z  +  a3 


b^Z  +  bgZ  +  bj 


is  constrained  by 


where 


i)  C.Z^  +  C,Z  +  C*  has  two  roots  inside?  the  unit 
1  L  i  circle. 

4  3  2 

1i)  D,Z  +  D?Z  +  D,Z  +D?Z  +  D,  has  two  roots  inside 
£  ~  '  the  unit  circle 


C1  =  a«  +  bl 
C2  =  a2  +  b2 

C3  =  a3  +  b3 

D-|  =(1/2)  (a3b1+a]b3) 

02  =(l/2)(a} b2+a2b3+a2b1  +  a3b2) 
°3  =  albl  +  a2b3  +  a3b3 


The  Adaptive  Control  (Tuning)  Mechanism  The  objective  of  the 
tuning  mechanism  is  to  have  the  plant  behave  like  the  reference  model.  It 
is  clear,  however,  because  the  positivity  constraints  preserve  the  GH^ 
loop  stability,  that  one  is  free  to  choose  the  tuning  law  as  long  as  it 
generates  bounded  gains. 


The  general  form  of  the  tuning  mechanism  resembles  the  generalized- 
steepest  descent  rule: 

Change  in  gain  =  -  X’  S  ’ 

where 

A 

X=  learning  gain  (possibly  time-varying  but  always  greater  than 
zero) 

A 

$=  positive  definite  normalization  factor  (possibly  time  varying). 

.  A 

gradient  related  direction 

One  may  choose  these  variables,  for  example,  based  on  least  squares  [15], 
gradient  methods  [16],  or  correlation/correlation-with-delay  [17].  We 
prefer  the  cress  correlation  law  because  it  has  previously  been  shown  to 
have  some  very  desirable  convergence  [18]  and  stability  properties  [19]? 
it  can  be  shown  to  generate  a  bounded  sequence  when  a  positive  real  refer¬ 
ence  model  is  used*  [22];  and  it  is  recursive  and  Is  easily  computed  on-line 
in  real  ti~ie. 

The  signal  in  the  cross-correlation  law  is  computed  as 
[  [I,  ©eT]  -g'  [I2  ©S]  dT 

where: 

a  =  [a..]  is  mx2  matrix  of  adaptive  gains. 

*  J 

e  is  the  plant-model  error  vector,  rxl 

S  is  the  vector  of  signals  entering/controll ed-by.  adaptive  gains 

*  We  can  show  this  even  for  a  mismatch  between  the  order  of  the  reference 
model  and  the  plant  and  with  insufficient  control  degrees  of  freedom.  Also 
we  conjecture  that  the  requirement  for  a  positive  real  reference  model  can  be 
eliminated  via  embedding. 


8  cc  «  . ,  2  2  i  \ 

3~2  is  the  m  x  i  matrix  of  m  x  i  partions  ) 

raij( 

I-j  is  the  identity  matrix  of  order  r,  the  length  of  e 
I2  is  the  identity  matrix  of  order  m,  the  column  dimension  of  <Z . 
®is  the  outer  (or  Kronecker)  product. 


The  outer  product  of  a  txm  matrix  M  with  a  pxq  matrix  N  is  defined 

as  the  x  mq  matrix  consisting  of  *xm  partitions  m..N  . 

^  J 


rw 

Ml 

2-7 

2-n] 

l  3*7 

_  i'll" 

■  5*7  ' 

'5'iT] 
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11 

u  14 
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Note  that  in  general  the  outer  product  generates  a  matrix  of  matrices. 

This  occurs  because  one  is  generally  minimizing  a  vector  of  error  signals 
and  there  is  a  matrix  of  parameters  with  which  one  could  conceivably  adapt 
(one  parameter  vector  for  each  element  of  the  transfer  function  matrix). 
Considerable  simplification  results,  however,  if  one  permits  no  coupling 
and  a  priori  selects  the  order  of  the  numerator  and  denominator  dynamics  of 
each  element  of  the  transfer  matrix.  The  resulting  GR  in  this  case,  extremely 
sparse. 


Parameterization  Of  The  Adaptive  Compensators  Figure  4.3  has  a 
more  detailed  block  diagram  of  the  adaptive  design.  The  precompensator  and 
feedback  compensator  are  factored  into  two  operators,  the  basis  transfer 
matrix  and  the  gain  matrix.  The  first  contains  all  the  aM  vat  i  ve/or  delay 
dynamics.  The  latter  contains  only  the  adaptive  gains.  The  realization  is 
possible  as  any  linear  finite  order  differential  (delay)  operator  H  that 
is  controllable  can  be  put  in  the  form  [20] 

H  =  N  P’1 
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Figure  4.3.  Positivity  Constrained  MRAC  Configuration 


Where  N  and  P  are  matrices  whose  elements  are  polynomial  functions  of  the 
derivative  (delay)  operator  D.  Our  adaptive  controller  will  vary  only  the 
compensator  zeros;  hence,  N  contains  all  the  adaptive  gains  and  P  is  fixed. 

N  can  be  now  factored  into  the  form  &>'  N  where  zt  contains  only  the  adaptive 

A 

/i  a  l  3  n/i  M  /•  3  >  r.nl  is  t  imo  i  nw3v*i  anf  Anovaf  a»*c 

«  I  I  •/  U  HU  It  -V  W  •  *  VIA  *  »  •  V  *  I  I  J  W  KIIV  a  II  •  Ul  I  UM  V  VI  U  W>  «*• 


Exampl e: 


*  -1 

Hence  if  N  and  P  are  combined  into  one  operator  (what  we  call  the  basis 
transfer  matrix)  all  the  adaptive  gains  are  left  in  one  gain  matrix.  Note 
that  if  H  is  proper,  derivative  operations  need  not  be  performed. 

On-Line  Computational  Requirements  Despite  the  complex  looking 
system  topology  shown  in  Figure  4-3,  the  computational  requirements  for 
digital  LSS  control  applications  are  modest.  Specifically,  if  one  allocates 
one  microprocessor  to  each  physical  control  station,  it  can  be  shown  that 

1)  The  decentralized  nth  order  digital  reference  modes  requires 

2n  -  multiplies 

2n  -  adds 

4n+2  -  storage  locations 
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2)  The  adaptive  law  with  the  basis  in  H0  equal  to  the  identity, 
adaptive  diagonal  forward  gains,  and  adaptive  diagonal  rate 
output  gains  requires 

9  -  multiplies 

1  -  divide 

25  -  storage  locations 

3)  The  adaptive  law  for  r  2nc*-order  digital  filters  (either 
in  H.|  or  H^)  requires 

1  +  11 r  -  multiples 
1  -  divide 

1  +  7r  -  adds 

7  +  1 9r  -  storage  locations 

4)  The  positivity  test  for  each  2nc*  *  order  digital  filter 
section  requires 

30  -  multiplies 

3  -  divides 

18  -  adds 

73  -  storage  locations 

EXAMPIE:  Assume  the  reference  model  is  6th  order,  the  forward  gains  and 
rate  feedback  gains  are  adaptive,  the  dynamic  feedback  is  an  adaptive  6th 
order  filter  (three  second  order  filter),  and  the  maximum  acceptable  compu¬ 
tational  time  delay  is  3  msec.  The  computational  requirements  per  control 
cycle  are: 

145  -  Multiplies 
11  -  Di-ides 

95  -  Adds 

664  -  Storage  locations  (safety  factor  of  two) 


If  one  further  assumes  that  (timewise) 


1  multiply  operation  =  3  add  operations 


Anri 


1  divide  operation  ~  3  multiply  operations,  then  the 
microprocessor  must  be  capable  of  performing  multiplication  in  about 
14  ^sec.  This  is  well  within  the  present  state-or'-the-art. 


4.2  DESIGN  COMPARISONS  USING  THE  DRAPER  TETRAHEDRON 


This  section  will  use  the  Draper  Tetraheurai  Truss  to  compare  fixed 
gain  controls  with  direct  adapti.e  controls  and  with  identification  followed 
by  fixed  gain  controls  (indirect  adaptive  control).  First  the  example 
structure  will  be  described;  then,  design  and  simulation  results  for  each 
approach  will  be  summarized;  finally,  a  comparison  of  computational  and 
hardware  requirements  will  be  made. 


4.2.1  Description  of  Draper  Tetrahedron 

The  Tetrahedral  Truss,  shown  in  Figure  4-4,  has  been  selected  as  a 
test  example  for  comparing  the  various  design  approaches.  This  structure 
was  devised  by  Draper  Labs  [21]  as  one  of  the  simplest  non-planar  geometries 
capable  of  representing  a  large  space  structure.  Despite  its  apparent 
simplicity,  it  models  the  feed-tower  in  a  generic  class  of  large  antenna 
applications.  The  tetrahedral  apex  represents  the  antenna  feed,  members  1-6 
are  part  of  the  feed  support  structure,  and  bi-peds  7-8,  9-10  and  11-12  are 
supports/controls  attached  to  an  inertially  stabilized  (assumed)  antenna 
dish.  Because  of  the  nature  of  the  real  problem,  the  motion  of  the  feed 
(apex)  must  be  controlled  through  actuators  and  sensors  placed  at  the  bi-peds 
none  can  be  placed  at  the  feed.  The  latter  constitutes  the  control  problem  ; 
i.e.,  control  the  feed  precisely  despite  large  modeling  uncertainty  and 
without  resorting  to  direct  measurement  or  control  of  its  motion. 
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The  fin-!  element  method  was  used  to  model  the  tetrahedron. 
Referring  to  Figure  4*4,  the  model  contains  ten  nodes,  each  with  three 
degrees  of  freedom,  and  twelve  truss  members.  Node  coordinates  and  element 
connectiveness  for  the  "nominal"  model  ate  listed  in  Table  4.1.  These 
elements  are  assumed  to  be  capable  of  resisting  axial  forces  only.  The 
masses  were  assumed  to  be  lumped  at  the  nodes  one  through  four. 

Six  co~located  sensors/actuators  were  assumed  to  act  as  member 
dampers  in  parallel  with  members  7  through  12.  The  member  dampers  were 
assumed  capable  of  sensing  relative  position  and  velocity  and  exerting 
force  in  the  axial  direction  only. 

The  objective  was  to  design  a  robust  controller  and  evaluate  the 
resulting  controller  on  the  full  order  nominal  model.  A  second  equation 
was  then  conducted  using  the  true  tetrahedral  system  parameters  s.<jwn  in 
Table  4.2.  Both  evaluations  were  based  on  the  intial  conditions  given  in 
Table  4.3.  The  goal  was  to  damp  the  1 ine-of-sight  motion  of  Node  1  in  the 
x  and  y  directions  to  less  than  .0004  and  .00025  units,  respectively,  in 
20  seconds. 


4.2.2  Digital  Fixed  Gain  Control  Of  The  Tetrahedron 


A  fixed  gain  control  system  for  the  tetrahedron  will  be  designed 
first,  to  serve  as  representative  of  a  standard  control  system  design.  The 
fixed  gain  control  system  design  for  the  tetrahedron  was  conducted  using 
the  general  frequency  domain/positivity  design  procedure  outlined  in  sub¬ 
section  2.1.6.  The  exception  was  to  use  the  digital  controller  design 
modification,  outlined  in  subsection  4.1.2,  so  that  comparisons  with  the 
digital  adaptive  design  are  fair. 


The  sampling  rate  was  chosen  to  be  .2  seconds.  This  rate  is  approxi¬ 
mately  tnree  times  faster  than  the  rate  required  to  reconstruct  the  highest 
controlled  mode  (.909  Hz).  At  this  sample  rate  one  can  compute  the  positi¬ 
vity  index  of  the  plant,  in  either  the  Z-domain: 


{.)  = 


•  4 

nn  I  2 


[G(2) 


G*(z  )] 


j-  T 


z=e 


sample 


[0 ,2*r ) 
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Table  4.1  Tetrahedron  Structural  Model  Parameters  (Nominal) 


Node  Point  Locations 


Node  X 


Y 


1 

2 

3 

4 

5 
5 

7 

8 

9 

10 


0.0 

0.0 

-5.0 

-2.887 

5.0 

-2.887 

0.0 

5.7735 

-6.0 

-1.1547 

-4.0 

-4.6188 

4.0 

-4.6188 

6.0 

-1 .1547 

2.0 

5.7735 

-2.0 

5.7735 

Restraints  X,Y,  4  Z  at  nodes  5-10 


Z 

10.165 

2.0 

2.0 

2.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


Element  Connectivities  and  Properties 


Material 


Element 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
n 
12 


Node  1 

1 

1 

1 

2 

3 

2 

2 

2 

3 

3 

4 
4 


Node  2 

2 

3 

4 

3 

4 

4 

5 

6 

7 

8 

9 

10 


Area 

1000. 

100. 

100. 

1000. 

1000. 

1000. 

100. 

100. 

100. 

100. 

100. 

100. 


Lumped  Masses 


Node  Mass 


1 

2 

3 

4 


2.0 

2.C 

2.0 

2.0 
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Table  4.2  Tetrahedron  Parameter  Variation 


Masses 


Node  Mass 

1  4.0 


Element  P»*operty  Changes 


Element 

Area 

1 

1200.0 

2 

150.0 

3 

150.0 

4 

'200.0 

5 

1200.0 

5 

1200.0 

7 

150.0 

8 

150.0 

9 

150.0 

10 

150.0 

11 

150.0 

12 

150.0 

Table  4.3  Initial  Conditions  Fc  -ontrol  Design  Evaluations 


Mode 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 


Displacement  (n) 


-.001 

-.003 

.006 

.01 

.001 

.03 

-.009 

-.02 

.008 

.02 

-.001 

-.02 

-.002 

-.003 

.002 

.004 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 
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or  in  the  W-domain*: 


6{a))  =  min  {  f[G(W)  +  j  ’  Wc  tO,  j-) 

The  characteristic  gains  required  in  the  frequency  domain  part  of 
the  design  can  similarly  be  computed  in  either  the  2  or  W  domain. 

The  decentralized  controller  resulting  from  a  design  in  the 
W- domain  was 

u  I  .Til  492  ♦  3,9211622  -  3.9863812  1 
H  6*6  l  ?  -  , ,  7092 lz  *  .8465, 8J  (CM) 

Figure  4.5a-f  show  the  transient  responses  of  the  nominal  and  per¬ 
turbed  systems  to  the  initial  conditions  given  in  Table  4.3  Figures  4.5a 
and  d  show  that  the  residual  modes  are  stable  despite  the  use  of  the  low 
order  decentralized  controller  (CN).  Figure  4.5b  and  e  show  that  the 
nominal  system  meets  design  specification  but  that  there  is  a  small  loss 
of  performance  using  the  perturbed  model.  Figures  4.5c  and  f  show  tha4- 
actuator  command  forces  are  modest. 

4.2.3  Digital  Direct  Adaptive  Control  Of  The  Tetrahedron.  The  next  step 
is  to  design  a  direct  adaptive  controller  so  that  one  can  compare  this 
approach  with  the  fixed-gain  approach.  The  digital  direct  adaptive  controller 
for  the  Tetrahedron  was  designed  using  the  positivity  constrained  adaptive 
control  theory  outlined  in  Section  4.1.  A  design  of  this  type  is  semi¬ 
automatic  in  the  sense  that  once  one  chooses  a  reference  model  and  controller 
configuration,  the  system  seeks  optimum**  performance  through  self-adjustment. 
The  design  effort  for  the  Tetrahedron  was  therefore  spent  on  choosing  a 
controller  configuration  which  could  be  compared  fairly  with  the  fixed-gain 

*  G(W)  =  G ( z  =  ) 

**  If  the  plant  were  of  finite  order  and  the  controller  had  sufficient  degrees 
of  freedom  to  make  the  model-plant  error  identically  zero,  the  optimum 
reached  is  the  same  as  in  an  absolutely  cooperative  game.  Since  our  situa¬ 
tion  violates  both  these  conditions,  the  game  theoretic  equilibrium  the 
controller  reaches  must  be  investigated  in  future  research. 
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Figure  4.5  Transient  Response  0f  Tetrahedron  Using  A  Digital  Fixed  Gain  Controller 


Refering  to  Figure  4.6,  the  sample  rate  was  chosen  to  be  the  same 
as  was  used  in  the  fixc'd  gain  controller  design,  .2  seconds.  The  reference 
model  picked  uses  the  eight  mode  (number  of  controlled  modes), 
nom i na 1  Tetrahedron  model,  but,  with  its  damping  augmented  to  reflect  the 
effects  of  the  nominal  fixed-gain  controller  (NC)  The  adaptive  decentral¬ 
ized  feedback  compensator  was  chosen  to  have  exactly  the  same  form  as  the 
controller  (NC);  i.e.,  has  the  same  poles  and  only  the  numerator  dynamics 
have  been  made  adaptive.  The  adaptive  precompensator  consists  of  a  diagonal 
gain  matrix,  to  compensate  for  "averaged"  mode  shape  uncertainty,  and  is 
only  required  during  "servo"  control. 


TH 

Finally,  the  cross-correlation  adaptive  law  used  in  the  i  path 
(each  of  six  control  stations)  is  computed  as 

-  02nT  ei^n>  *  [Signal  entering  k.(n)] 
k.(n+l)  =  k. (n)  -  10  e  1  _2 _ 1 . 

1  1  f-  r 


a^n+1)  -  a.(n)  -  10  e 


-.02nT  ei^n>  *  CSign?' 


V'n)l 


1  +  NO. 


-  o?nT  Mn)*  [Signal  entering  b.-(n)] 

Mn+1)  =  b.(n)  -  10  e  -^ni  _J _ 1 

1  +  NORM. 


where 

NORM.(n)  *  [Signal  Entering  k^(n}]^  +  [Signal  Entering  a^(n)]^ 

2 

+  [Signal  Entering  b . ( n ) 3 

e.(n)  c  Plant  output^n)  -  Desired  output. (n) 

Note,  the  exponentially  decaying  learning  gains  are  used  to  reduce  the 
effects  of  the  non-matchable  parts  of  the  controlled  plant  and  reference  model. 

Figure  4,u  shows  a  block  diagram  of  the  adaptive  controller  with  the 
full  order  (12  mode)  perturbed,  continuous  dynamics  of  the  Tetrahedron.  The  node 
of  operation  is  to  have  regularly  scheduled  "tuning"  periods  during  which 
a  training  command  is  injected  into  the  system  and  the  adaptation  loop  is 
turned  cn.  Normal  operation  follows  with  the  adaptive  loop  turned  off  and 
the  control  gains  remaining  fixed  at  their  last,  converged  value.  This  is 
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important  from  the  viewpoint  that  the  training  signals  required  for  adapta¬ 
tion  can  be  annoying  during  a  regulation  mode.  The  proposed  mode  of  operation 
would  turn  this  signal  on  and  off  together  with  the  adaptive  loop. 

Figures  t.7a-d  show  the  time  history  of  the  various  adaptive  system 
signals  during  a  "tuning"  session,  as  simulated  on  a  digital  computer. 

Figure  4.7a  shows  the  training  signals  input  to  stations  1  and  2  of  leg  1. 

It  is  a  periodic  sequence  consisting  of  5  second  pulses  followed  by  a  rest 
period  of  17.5  seconds.  The  training  system  signals  at  the  other  legs  were 
of  similar  periodicity  but  phased  so  that  the  apex  was  gyrated.  Figure  4.7b 
shows  the  values/change  in  values*  for  the  adaptive  gains  during  this  same 
period.  The  initial  values  for  these  gains  were  assumed  to  be  at  the 
nominal  controller  values  (NC).  Figures  4,7c  and  d  compare  the  reference 
and  "true"  response  of  modes  1  and  5  during  the  training  session.  Mode  5 
is  shown  to  improve  slightly  but  the  effect  on  mode  1  is  imperceptable.  It 
can  be  shown  that  the  reason  for  this  poor  state  tracking  of  mode  1  is  due 
to  its  poor  observability/controllabil ity  from  the  tetrahedral  legs.  Indeed, 
analysis  reveals  that  modes  1  and  2  contribute  25-50  db  less  nower  than 
other  mo.es  observable  at  the  legs.  Hence  these  modes  remain  "in-lhe-noise" 
with  rther  unmodelled  residual  modes.  This  effect  makes  it  difficult  for 
the  adaptation  algorithm  to  extract,  any  usef  il  information  from  the  output 
error  signal.  The  net  result  presents  a  problem  as  the  apex  motion  is 
dominated  by  these  poorly  observable  modes. 

Figure  4.8  shows  the  transient  response  of  the  perturbed  system  with 
the  adaptive  gain  held  fixed  at  the  converged,  post  tunning  session  values. 
The  adaptive  control  law  performed  well.  But  due  to  the  Fact  that  the  modes 
that  affect  motion  at  tne  tetrahedron  apex  are  poorly  observable  from  the 
sensor  measurements  in  the  bi-ped  legs,  the  adaptive  law  did  not  perform 
any  better  than  the  (very  robust)  fixed-gain  law  that  it  is  being  compared 
with.  It  should  be  noted,  that  the  riur.-aaaptive  law,  due  to  its  robustness, 
was  able  to  just  about  meet  the  20  second  settling  specifications  even  when 
the  perturbed  structural  modal  was  used. 


*  Forward  and  rate  feedback  gains  are  measured  relative  to  zero.  Filter 
gains  are  measured  relative  to  nominal  filter  gain  values. 
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Comparing  these  results  with  those  from  Section  4.2,  it  can  be  con¬ 
cluded  that  adaptive  control  is  feasible  and  will  work  well  for  LSS,  but 
that  full  advantage  of  its  power  can  only  be  taken  in  situations  where  (i)  the 
motion  at  the  point(s)  of  interest  on  the  structure  is  measurable  (or  can 
be  estimated)  in  real  time,  and  (ii)  changes  in  the  dynamic  model  of  the 
structure  are  sufficiently  drastic  and  occur  in  real  time  (e.g.,  on-orbit 
reorientation  of  articulated  substructures) ,  which  would  cause  problems  for 
fixed-gain  controllers. 

4.2.4  Identification  of  the  Tetrahedron  Followed  By  Fixed-Gain  Digital 
Control  [22]. 

After  having  designed  and  compared  fixed-gain  and  adaptive  controls, 
the  third  alternative,  identification  followed  by  fixed-gain  control,  will 
be  designed  and  compared.  The  maximum  likelihood  identification  approach 
described  in  Section  3  can  be  combined  with  the  fixed-gain  control  technique 
outlined  in  Subsection  2.1.6  and  4.2.2  to  design  an  indirect  adaptive 
controller  for  the  Tetrahedron.  In  the  first  part,  system  identification,  a 
physically  parameterized  finite-dimensional  state-space  model  is  derived 
through  the  use  of  the  finite  element  method.  The  parameters  to  be  estimated 
are  the  cross-sectional  area  of  all  12  tetrahedral  truss  members  and  the 
weight  of  the  masses  located  at  nodes  1  through  4.  There  are,  therefore,  16 
parameters  to  be  estimated.  Note  that  if  modal  frequencies  and  the  corres¬ 
ponding  mode  shapes  had  been  chosen  to  parameterize  the  system,  the  number  of 
unknown  parameters  would  be  56  for  an  8  mode  model,  and  84  for  a  12  mode 
model.  Hence,  by  identifying  physical  parameters,  one  has,  in  fact,  reduced 
the  number  of  unknown  parameters. 

Figure  4.9  shows  the  input  signals  used  to  excite  the  structure  and 
generate  the  data  required  by  the  maximum  likelihood  algorithm.  The  signals 
are  pulses  of  magnitude  1  and  last  .25  seconds.  Each  pulse  is  staggered  in 
time  and  is  applied  to  the  actuators  located  at  the  truss  members  identified 
in  the  figure.  The  rate  sensor  responses  to  these  pulses,  which  are  referred 
to  as  the  Tetrahedral  data,  were  simulated  using  the  nominal  model*,  and, 
random  noise  was  added  (  a  =  .002).  The  perturbed"  model  parameters  were 
then  chosen  as  the  initial  parameter  values  used  in  the  identification 

*  The  role  of  the  nominal  and  perturbed  models  were  inadvertently  reversed 
during  the  identification  procedure.  However,  this  does  not  impact  the 
results  to  be  presented. 
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Figure  4.9  Input  Signals  Used  To  Exite  Structure  During 
Identification  Procedure. 
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algorithm. 

Table  4.4  shews  the  convergence  of  the  parameter  values  In  the 
Newton-Raphson  iteration,  using  5  seconds  of  -  All  15  parameters  con¬ 
verged  in  9  iterations  and  the  final  estimates  are  very  close  to  the  true 
parameter  values. 

Table  4.5  shows  the  equivalent  modal  quantities  represented  by  the 
physical  parameters  in  Table  4.4.  These  are  obtained  by  finite  element 
processing  of  each  physical  parameter  set  and  demonstrate  that  indirect 
(i.e.,  via  physical  parameter  estimation)  modal  parameter  estimation  is 
al so  excellent. 


Figure  4.10  compares  the  time  histories  of  the  Tetrahedron  data  and 
the  estimated  data  as  obtained  by  using  the  final  parameter  estimates  given 
in  Table  4.4.  It  can  be  seen  the  fit  is  very  good  and  it  can  be  shown  that 
the  error  consists  primarily  of  measurement  noise. 

The  second  part  of  the  indirect  adaptive  approach  Is  control  system 
design  tune-ups.  The  feedback  law  tune-up  was  based  on  the  W-domain  fixed- 
gain  design  approach  described  earlier.  In  order  to  make  comparisons  with 
the  direct  adaptive  design,  however,  the  poles  of  the  updated  controller 
were  fixed  at  the  nominal  controller  values  and  only  the  numerator  of  the 
controller  was  updated  (i.e.,  recomputed).* 


The  feedback  law  based  on  the  final  parameter  estimates  is: 


..  ,  .  6-9611 6z  -  7.07695 

H=  IKvfi  .  14.3567  r -x - 

bxb  z  -1.70921 z+. 846578 


Note  that  only  the  numerator  of  (NC)  was  modified  to  create  (NU). 

Figure  4.17  shows  the  transient  response  of  the  perturbed  Tetrahedron 
to  the  initial  conditions  given  in  Table  4.3.  The  traces  show  that  all 
specifications  are  met  by  using  the  updated  controller. 

*  Adaptive  design  theory  permits  adaptive  "pole  motion",  but  not  simul¬ 
taneously  with  "zero  motion". 
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Table  4-4  Convergence  of  Physical  Parameters  In  The  Newton-Raphson 
Iteration. 


Initial 

Parameters  Guess 


Iterations 


Final 

Estimate 


9 


ai 

1200 

1177 

1290 

1230 

1067 

1005 

1000 

a2 

150 

149 

190 

67 

93 

101 

100 

a  *  l 
i  1 

1  CP 

167 

75 

53 

107 

103 

100 

ra 

(V 

a4 

1200 

1257 

1103 

1006 

1017 

1008 

1000 

L 

< 

a5 

1200 

1056 

1078 

1135 

1041 

1030 

1000 

c 

o 

•e* 

a6 

1200 

1168 

1138 

1068 

1061 

999 

1000 

u 

a; 

CO 

i 

a7 

150 

105 

104 

96 

105 

101 

100 

</> 

O 

a8 

150 

190 

296 

188 

105 

101 

100 

O 

1 

l- 

<U 

-o 

a9 

150 

145 

132 

111 

99 

98 

100 

3: 

aio 

150 

150 

145 

129 

101 

100 

100 

</> 

</> 

all 

150 

134 

112 

103 

93 

94 

100 

— j 

V. 

»— 

al  2 

150 

15* 

159 

107 

98 

98 

100 

l/m] 

.250 

.256 

.293 

.364 

.441 

.483 

l/m2 

.500 

.476 

.438 

.453 

.462 

.483 

%A 

</> 

l/m3 

.500 

.481 

.420 

.441 

.508 

.500 

z: 

l/m4 

.500 

.474 

.444 

.435 

.511 

.515 

Table  4-5  Convergence  of  Equivalent  Modal  Parameters  During  The 
Newton-Raphson  Iteration. 


MODE 

FREQUENCIES* 

(Hz) 

INITIAL 

ITERATIONS 

FINAL 

True 

0 

1 

3 

5 

7 

9 

“1 

.186 

.190 

.168 

.154 

.201 

.213 

.214 

“2 

.233 

.235 

.241 

.206 

.251 

.263 

.265 

w3 

.471 

.459 

.460 

.416 

.447 

.457 

.460 

“4 

.566 

.533 

.495 

.475 

.473 

.471 

.471 

u5 

.612 

.602 

.598 

.567 

.530 

.537 

.541 

“6 

.819 

.798 

.784 

.712 

.666 

.665 

.669 

w7 

.903 

.882 

.818 

.721 

.743 

.743 

.742 

“8 

.908 

.887 

.839 

.745 

.755  ■ 

.758 

.757 

U9 

1.422 

1 .378 

1 .366 

1 .360 

1 .352 

1 .359 

1.359 

w10 

1.639 

1 .609 

1 .510 

1 .458 

1.472 

1 .472 

1.472 

“11 

1.738 

1 .639 

1 .639 

1  .674 

1.649 

1  .649 

1 .637 

u12 

2.222 

2.155 

2.039 

2.044 

2.050 

2.054 

2.054 

*  Mode  frequencies  were  computed  from  estimated  physical  parameters. 
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Table  4.6  compares  the  computational  and  hardware  requirements  for 
the  Tetrahedron.  The  table  shows  that  the  steady-state,  on-orbit  control 
computation  requirements  is  the  same  for  the  three  approaches  studied.  The 
difference  occurs  only  during  the  training  or  learning  session  required  when 
using  the  direct  adaptive  approach.  During  this  period,  the  computational 
requirements  increase  to  account  for  the  adaptive  control  law  and  the  on-line 
positivity  tests  that  ensure  stability.  In  terms  of  required  on-orbit  hard¬ 
ware,  fixed-gain  controllers  are  the  simplest,  and,  depending  on  the 
acceptable  computational  delays,  can  be  implemented  using  one  microprocessor. 
Similarly,  decentralized  fixed-gain  control  may  require  only  one  microprocessor, 
but,  a  storage  device  must  also  be  provided  to  store  the  input-output  data 
required  for  the  parameter  identification  processing  to  be  done  on  the 
ground.  The  direct  adaptive  controller  would  require  six  microprocessors 
unless  processors  with  1-2  psec  multiply  instruction  cycles  are  used. 


4.3  Conclusions  and  Recommendations 

Chapter  4  has  compared  adaptive  and  non-adaptive  controls.  It  has 
been  shown  that  positivity  theory  can  be  used  to  ensure  that  either  type  is 
stable  and  that  each  type  is  realizable  with  current  hardware.  Fixed  gain 
controllers  were  found  to  yield  sub-optimal  performance.  It  was  uncovered, 
however,  that  direct  adaptive  controls  do  not  significantly  improve  the 
fixed  gain  control  performance  if  the  controlled  quantities  are  weakly 
observable/controllable  from  the  "available"  sensors  and  actuators.  Indirect 
adaptive  control,  where  one  identifies  physical  parameters,  can  circumvent 
this  difficulty  if  one  accounts  for  the  weak  observability/ control  lability 
in  the  design  phase  and  uses  any  other  available  information. 

It  is  recommended  that  further  studies  De  conducted  in  the  following 

areas : 


Determine  if  the  adaptive  control  reduced  the  bandwidths 
requirements  over  fixed-gain  controllers. 


Table  4*6  Comparison  of  Computational  and  Hardware  Requirements  For  The  Tetrahedron 
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FOR  EACH  OF  6  CONTROL  CHANNELS 


Explore  control  design  tradeoffs  when  LSS  dynamics  are 
time-varying  very  rapidly  with  respect-to  the  time  used 
for  parameter  estimation  and  controller  updating. 

Determine  if  weakly  observable  modes  that  must  be  controlled 
can  be  identified  using  modal  parameter!' zations  rather  than 
physical  parameter izat ions. 

Design  adaptation  algorithms  which  "optimize"  the  basis- 
transfer-matrix  elements  on-line  (see  Figure  4.3). 
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APPENDIX  A.  Draper  Example  Structure  Model  #2 


The  Draper  Model  § 2  shown  in  Figure  A-l  is  a  large  space  structure 
model  provided  by  R.  Strunce  and  T.  Henderson  of  Charles  Stark  Draper 
Laboratory  (CSDL)  for  the  purpose  of  testing  and  eva^ating  controller 
design  methodologies  and  system  identification  algorithms.  The  structure 
consists  of  a  flexible  optical  support  structure  ("clean  structure")  and 
an  equipment  section  ("dirty  box" )y  separated  from  the  "clean  structure" 
by  passive  isolators.  The  optical  support  consists  of  the  primary  and 
tertiary  mirror  r"pport  truss,  and  a  lower  support  truss  for  the  sec¬ 
ondary  mirror  and  focal  plane  sensor.  The  line-of-sight'(LOS)  path  is 
also  shown  in  Figure  A-l.  The  isolated  equipment  section  consists  of  a 
central  rigid  body  containing  the  spacecraft  attitude  control  system. 

The  equipment  section  also  supports  two  symmetric  solar  panels. 

Draper  has  provided  one  nominal  model  and  two  perturbed  models 
(nominal  model  witn  parameter  variations),  and  an  expression  describing 
the  LOS  error,  presented  in  Section  A. 3  of  this  Appendix. 

A.l .  Nominal  Model 

The  finite  element  model  is  described  in  Figures  A-2  through  5  and 
Table  A-l.  The  modal  frequencies  and  mode  shapes  are  listed  in  Table  A-2. 
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Figure  A-l .  Draper  Example  Structure  Model  #2 


X?  i'\ 
Xv 


\  \ 


^x.  X 

'  V 


NODE 

X 

Y 

26 

-5.0 

o 

o 

27 

-4.0 

3.0 

20 

-4.0 

-3.0 

29 

4.0 

3.0 

30 

4.0 

-3.0 

31 

5.0 

o 

© 

52 

-4.0 

10.0 

33 

4.0 

10.0 

34 

-4.0 

-10.0 

35 

4.0 

-10.0 

36 

-4.0 

3.0 

37 

-4.0 

o 

1 

)B 

4.0 

3.0 

39 

4  0 

-3.0 

DIAGONAL  MEMBERS  OMITTED  TOR  CLARITY 


MEMBER  • 


MEMBER  • 


Figure  A- :  Nominal  Model  -  Upper  Support  Truss 


WMiSTiiMtiV  L  MM 


Table  A-l.  Nominal  Model  Member  Properties 


TYPE 

200 

TYPE 

300 

20  cm 

dia.  x  .1  cm 

Round  Tube 

20 

cm 

dia  x 

.05  cm  Round  Tube 

A  = 

6.250  x  10*4 

m2 

A 

s 

3.133 

x  10*4  m2 

I  = 

3.075  x  10'6 

m4 

I 

a 

1.559 

x  10  6  m4 

J  » 

6.189  x  10~6 

4 

m 

J 

a 

3.118 

x  10  £  m4 

TYPE 

400 

TYPE 

500 

25  err 

i  dia  x  .05  cm 

Round  Tube 

40 

cm 

dia  x 

.075  cm  Round  Tube 

A  * 

3.919  x  10'4 

m2 

A 

r 

9.407 

x  10"4  m2 

I  * 

3.0496  x  10"£ 

1  m4 

I 

a 

1.874 

x  10'5  m4 

J  * 

6.099  x  10"6 

m4 

J 

a 

3.749 

x  10-5  m4 

Table  A-2.  Nominal  Model  Modal  Frequencies  and  Mode  Shapes 


Freq. 

(HZ) 

Description 

Freq. 

(Hz) 

Description 

1-6 

0.0 

rigid  body 

17 

1.72 

torsion 

7 

.  145 

isolator  Y-rotation 

18 

1.82 

2nd  S.P.  X-Z  plane 

8 

.263 

Isolator  Z-rotation 

19 

1.82 

2nd  S.P.  X-Y  plane 

9 

.317 

)st  S.P.  X-Z  plane 

20 

1.89 

1st  S.P.  torsion  asym 

10 

.  333 

1st  S.P.  X-Y  plane 

21 

2.36 

1st  bending 

11 

.443 

Isolator  Z-trans 

22 

2.99 

1st  S.P.  torsion  sym 

12 

.577 

Isolator  Y-trans 

23 

3.18 

3rd  S.P.  X-Y  plane 

13 

.  581 

Isolator  X-trans 

24 

3.39 

3rd  S.P.  X-Z  plane 

14 

1.22 

2nd  S.P.  X-7.  plan- 

25 

5.16 

2nd  S.P.  tors-  „n 

15 

1.  30 

2n<1  S.P.  X-Y  plan- 

26 

5.26 

2nd  S.P.  torsion 

16 

1.  35 

Isolator  X-rotation 

2i 

7.87 

3rd  S.P.  torsion 
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Table  A-2.  Nominal  Model  Modal  Frequencies  and  Mode  Shapes  (continued) 


A. 2.  Perturbed  Models 


Two  perturbed  models  are  given  -  Case  if 2  and  Case  #4.  The  changes 
in  masses  and  stiffnesses  are  listed  in  Table  A-3  and  A-4.  All  other 
parameters  are  unchanged. 


Table  A-3.  Changes  in  Masses  for  Perturbed  Models 


Node 

Original 

Mass 

Case  n 

Case  r 4 

71 

375.0 

350.0 

375.0 

23 

375.0 

350.0 

375.0 

29 

375.0 

350.0 

375.0 

30 

375.0 

350.0 

375.0 

32 

500.0 

500.0 

550.0 

33 

500.0 

550.0 

450.0 

34 

250.0 

300.0 

300.0 

35 

250.0 

250.0 

200.0 

Table  A-4.  Changes  in  Stiffnesses  for  Perturbed  Models 


Merab . 

NA 

NB 

Section  Property* 

Case  #2 

Case  £4 

76 

8 

14 

1 

1 

80 

9 

15 

1 

2 

81 

11 

17 

1 

1 

85 

13 

19 

1 

2 

99 

15 

32 

1 

2 

101 

17 

33 

1 

2 

+  Section  Properties 


*1 


-  25cm  DIA.  x  .10 
A  =  7.823  x 

I  =  6.063  x 

J  =  1.213  x 


cm  round  tube 
Kf4  m2 
lO"6  m4 

ur5  »4 


=2  -  25  cm  0IA.  x  .025  cm  round  tube 

A  =  1.962  xlO'4  m2 

i  =  1.529  x  10'6  m4 

0  =  3.059  x  10'6  m4 


A- 11 


A . 3 .  LOS  Error  Algorithm 

Equations  (A-2),  (A-3),  and  (A-4),  relate  optical  surface  motion 
and  the  resultant  LOS  rotations  about  the  X  and  Y  axes.  The  notation 
and  orientation  is  consistent  with  the  structural  analysis.  The  move¬ 
ment  of  each  mirror  is  defined  to  be  the  displacements  and  rotations  in 
the  global  X,  Y,  Z  directions  at  a  point  on  the  sphere  (primary,  tertiary) 
or  plane  (focal  plane)  which  intersects  the  optical  (global  Z)  axis. 

The  LOS  rotation  and  defocus  are  given  by 

LOSX  =  Y/F  LOSY  =  X/F  Defocus  =  Z 

where 

F  =  8.051 

X  =  A]  C-Xp  +  Xt  -  R-0Yp  +  A2*6Ys  -  2T-eYt]  +  Xt  -  Xf 

Y  =  M1  [-Y  +  Yt  +  R-eXp  -  A2.0Xs  +  2T.9Xt]  +  Yt  -  Yf 

Z  ■  A3  CZp  -  2Zs  *  V  +  Zt  -  Zf 

where 

A1  =  0.2S87 
A2  =  93.90 
A3  =  0.0892 
R  =  53.9 
T  »  66.95 

The  terms  X. ,  Y. ,  Z.  9X. ,  9Yj,  9Z.  refer  to  the  translations  and  rota¬ 
tions  in  the  global  X,  Y,  and  Z  directions  of  the  primary  (p)  and 
secondary  (s),  tertiary  (t),  and  focal  plane  (f). 

The  expresssions  for  the  translation  and  rotation  of  each  mirror 
in  global  coordinates  are  a  function  of  the  support  point  displacements. 

Primary: 

Xp  =  X34  +  1,25  Y34  '  1,25  Y35 
Yp-0*50  ^34  +  0.50  Y35 

Zp  =  -0.2143  Z34  +  -0.2143  Z35  +  0.7143  Z2g  +  0.7143  Z3Q 


(A-l) 

(A-2) 

(A-3) 

(A-4) 


A-l  2 


Secondary: 


6Xp  =  -0.0714  Z34  -  0.0714  Z35  +  0.0714  Z2g  +  0.0714  Z3Q 

eYn  =  0.0125  Z,,  -  0.125  Z,, 

P  34  35 

0Zn  =  -Q.125  Y,.  +  0.125  Z,c 

p  34  35 


o 

•tf* 

X 

II 

co 

X 

exs 

o 

**• 

X 

CD 

II 

Vs  *  Y40 

■T. 

=  9¥40 

Zs  ■  Z40 

'eZs 

9Z40 

Tertiary: 

Xt  =  X27  +  -°-3750  Y27  +  °*3750  Y29 
Yt  =  0.50  Y27  +  0.50  Y2g 

Zt  =  0.0714  Z2?  +  0.7143  Z2g  -  0.2143  Z32  -  0.2143  Z33 
0Xt  =  -0.0714  Z2?  -  0.0714  Z2g  +  0.0714  Z32  +  0.0714  Z33 
0Yt  =  0.125  Z2?  -  0.125  Z2g 
0Zt  =  -0.125  Y27  +  0.125  Y2g 

Focal  Plane: 

Xf  =  Xn  +  0.6250  Yn  -  -0.6250  Yg 
Yf  =  0.50  Yn  +  0.50  Yg 

Zf  =  Z40 

0Xf  =  0.10  Zn  +  0.10  Zq  -  0.20  Z4Q 
9Yf  =  -0.125  Zn  +  0.125  Zg 
0Zf  =  0.125  Y^  -  0.125  Yg 
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